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DISTRIBUTION OF POTENTIAL IN A TWO-LAYERED MEDIUM 

DUE TO AN INTERNAL SOURCE AND SINK AND THE DETER- 

MINATION OF THE APPROXIMATE AVERAGE RESISTIVITY 
OF THE MEDIUM! 


By LACHLAN GILCHRIST, NORMON ROSTOKER, AND BEN BERNHOLTZ 


Abstract 


A brief review is presented of simple forms of electrical layout which are used 
in resistivity investigations on the upper surface of a two-layered medium of 
which the crust of the earth is assumed to be an example. The associative formulas 
for resistivity are given. 

The necessity for equalization in magnitudes of the source and the sink in all 
experimental methods of resistivity investigation is emphasized in order to be 
in accord with the basis of development of the formula, and in order to eliminate 
experimentally pseudo anomalies which may appear in homogeneous media. 

Formulas are developed for use in conjunction with experimental results 
obtained by the use of drill holes whereby the source and sink are embedded in 
a two-layered medium. Distinction is made in the cases where the source and 
sink are remote from, or are in the neighborhood of, the plane of separation of 


the layers. 
Some experimental results on the values of the interbowl resistance, R, are 


presented in comparison with the values of R that would be obtained from a 
medium of two homogeneous layers of a suitably chosen resistivity p. The local 
anomalies in these layers appear clearly. 


Part I 
Introduction 


A. Extensive theoretical and experimental research on this problem has been 
carried out by numerous investigators when the Electrical Source and Sink 
have been placed on the upper surface of the uppermost layer of the medium. 
Much of the work presented carries the assumption of the existence of an 
electrical source alone. The accompanying electrical sink is supposed to be 
placed at a great distance from the source and the effect of the sink on the 
potentials at points near the source is assumed to be negligible. In several 
electrical layouts, such as the Wenner layout, the effect of both source and 
sink must be taken into account in the development of the formula for the 
determination of the apparent average specific resistance or resistivity, p, of 
the medium. The formulas which were developed initially for the deter- 
mination of p in these cases apply, in general, only to a homogeneous medium, 
and were modified subsequently for specified cases of inhomogeneity such as 
that of a horizontally or obliquely layered medium. 


1 Manuscript received in original form June 28, 1949, and, as revised, October 18, 1949. 


Contribution from the Department of Physics, University of Toronto, Toronto, Ont. 
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B. The initially developed formulas for the determination of p in a homo- 
geneous medium used with various surface experimental electrical layouts are 
as follows: 


1.a. The Schlumberger layout with two point current electrodes as a source 
C, and a sink C, 





2r 
— R a , (1) 
l 1 1 


— ee YH 


a b L-a L-b 
where R = measured interbow] resistance, 


a, b = distances of the source C; from the potential electrodes P;, Pe 
respectively, 
and LZ = distance between C;, and C:. 


If LZ —aand L — bare large then 


p= =... ’ (2) 
ee 

b 

which is the approximate formula applicable in the region adjacent to the 

source C; when C; is on the surface of a homogeneous medium. This electrical 

layout is shown in Fig. 1. 


ze 
a 


; : Zz 2L 
In Fig. 1, if a = and b = — the layout becomes that suggested by 


e 


Wenner in 1916 and the formula for the resistivity becomes p = 27Ra. 





1.6. The electrical layout with collinear source C, and two equal sinks C. and 
C, equidistant and on opposite sides of Ci: 
2r 
1 L 7 (3) 


1 
cin "ta, rntogeniins Stl: diacetate 
a 6 B-@ L-—# 





p=R 


This is a very useful modification of the Schlumberger layout. aandb 
have the significance given in Fig. 1. 


1l.c. The electrical layout with source C; and four symmetrically arranged equal 
sinks C2, C2’, C2!’, C2’. 


= R— - — - mc ie 
. 1 . oe es a (4) 
a 





2P-e 2 B-B 2 P+e 2 +e 





1 
b 
This layout was designed and has been used effectively to explore the region 


about a single point source C; as center. This arrangement is, however, not 
movable conveniently along a profile. 


| 





ee 
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Ammeter Generator 





C,, C2 - Source & Sink 
P, , P2 - Potential electrodes 
Fic. 1. Schlumberger layout. 


2. The line electrode layout which has a uniform central current line elec- 
trode C; as source and two symmetrically placed parallel uniform line elec- 
trodes C2 and C, as equal sinks. 

- : al ee ; 
For the region near Ci, p = R———, where / is the length of the line 
) 

log. — 
a 
electrode and R is the measured cylindrical interbowl resistance. This 
electrical layout is shown in Fig. 2. 


line electrode sink C2 






To Potentiometer 







Generator 


central line electrode source 
Ammeter 





line electrode sink C2 


Laren BE ood 


Fic. 2. Central line electrode layout. 


3. A modification of the central line current electrode layout that has been 
found useful in exploration in the field is shown in Fig. 3. 


is applicable for the region near the 





Tr e T. 
_ The formula for p, viz., p = - 


log, — 
a 


middle half of the line source C}. 


t 
{ 
: 
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4. It is essential in all electrical layouts that all sinks should be equal, and, 
further, that the source should be equal to the sink as applied to the actual 
region under exploration. A careful study of many investigations which have 
made use of these layouts reveals that this is usually overlooked or neglected. 


TRI 


, & 
loge 


© 
©) 


P, ni 
Ce L Ce 


C,-a line current electrode as source, length 2 
C,C;- point current electrodes as equal sinks 


$ - alternating current generator 
A- alternating current ammeter 
V- vacuum tube voltmeter or potentiometer 
P, P.- potential electrodes 
Fic. 3. Single line electrode layout. 


C. Since these exploratory layouts are used for the purpose of locating and 
delineating inhomogeneities or true anomalies in that part of the medium 
which is under exploration, it is evident that the divergences of the measured 
values of R, which are obtained experimentally by means of any of these 
electrical layouts, from the theoretically calculated values of R which are 
obtained from a suitably chosen value of p become thus of critical significance 
for interpretation. For example, in the case of an assumed layered medium, 
one or two strategic experimental values of the interbowl resistance, R, and 
of the resistivity, p, usually those adjacent to the source C; and therefore 
approximately applicable to shallow depths in the uppermost layer, should be 
made use of as the foundation for the subsequently calculated values of R 
and p for an extensive homogeneous medium. Interpretation of the remaining 
experimentally measured values of R and p becomes thus the interpretation 
of the divergences from the initially calculated values of R and p for a homo- 
geneous medium. 


It becomes obvious, however, that the theoretically calculated values of R 
from the initial p which have been developed on the basis of a homogeneous 
medium may be inadequate, and may render uncertain, or may at least make 
difficult, the comparison of the calculated and experimentally measured values 
of R and p even for a nonhomogeneous medium of a simple form. Several 
investigators have endeavored to obtain more precise formulas for various 
assumed forms of inhomogeneity in the medium. In particular, a great deal 
of attention has been given to the case of a medium with layered inhomo- 
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geneities, largely because the formulas which have been derived may be used 
with a degree of confidence in exploratory work on the earth since it is a type 
of geological structure which is approximately layered in many regions. The 
assumption is made in general that the source and sink are on the upper surface 
of the uppermost layer, and, further, many of the theoretical investigators 
have assumed the existence of a source only. 


D. If the medium is constituted of layers of material of resistivities p1, p2, then, 
the electrical source may be placed internally in the layers of the medium. 
This, for example, may be done in the case of the earth by the use of a drill hole 
into which a contact current electrode as a source C; may be lowered to various 
depths. 


Accompanying equal sinks C2, C2’, C2’’, C2'"’ are placed on the surface of the 
ground symmetrically about the source C, and at a considerable distance from 
it. The potential contact points P; and P, are also placed in the drill hole, 
either above or below, or both above and below the source Ci, and at such 
distance from C; as to avoid the steep gradient of potential adjacent to the 
contact surface of C;. The layout is shown in Fig. 4a. 


C2 bo north line electrode 
line electrode & line eiectrode C2 





C,,C2-Current electrodes P, , P. - Potential electrodes 
Collar is 1125 ft. below surface 


Fic. 4a. Set-up for electrical logging. 


If the equal sinks C2, C2’, C2’, C2’”’ are at a considerable distance from C; 
the distribution of potential and the resultant interbowl resistance R may be 
measured, and if the source C; is at considerable depth the initial resistivity 
may be obtained from the approximate formula for a homogeneous medium, 


viz. 
: An 


p=R = (5) 


b 


1 _ 
a 

As the source Ci, together with the potential electrodes, is raised in succes- 
sive steps near to the surface of the ground, i.e., electrical logging, the formula 
for resistivity must be modified and at the surface becomes for a homogeneous 
medium formula (4) section B with 4” replacing 27. It is apparent however 
that the formula p = 42 R ab/(b — a) may be used if C; is left at the bottom, 
and the potential electrodes are moved up the drill hole. 


a I eens 


5 NEA STARE EN RRA cn UTNE 


ec OR PREC E TON TITERS SNES ESE RS 
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If by means of switches the electrical current is sent from C, to the equal 
sinks C2, Cs’, C/’, Co’”’ in succession, it is apparent that the interbowl resis- 
tance R, and the resistivity p, in various directions about the drill hole may 
be compared. If the electrical current is sent symmetrically in all directions 
about the drill hole it is apparent that the drill hole undergoes a form of elec- 
trical logging and even materials at some distance from the drill hole will have 
an influence on the results of the “‘logging’’ for R and p. If these results for 
Rand pare compared with the results obtained by the well known Schlumberger 
method of electrical “‘logging’’, significant information on R and p for ma- 
terials remote from the drill hole may be obtained. 


FE. In the case of a layered medium both electrical source and electrical sink 
may be placed internally in the layers of the medium. This, for example, may 
be done in the case of the earth by the use of two neighboring drill holes into 
which are lowered both source and sink respectively, as shown in Fig. 5. 


Since these drill holes are in many cases close together it is generally not 
justifiable to neglect the effect of either the source or sink in the delineation 
of the distribution of potential in the medium, or to make use of the formula 


Ar : : hai 
p = R———— fora homogeneous medium to obtain the calculated initial values 
] 


a b 

of p and the subsequently calculated values of R which are used for purposes 
of comparison with the measured values of R. This must be replaced by the 
formula which was developed for a homogeneous medium in which the effect 
of both source and sink are taken into account, viz., formula (1) Section B 
with 4m replacing 2x. These formulas are of course applicable only if both 
source and sink are at a considerable depth below the upper layer. However, 
for a layered or other inhomogeneous medium this formula is inadequate, and 
the uncertainty previously indicated in the comparison of the calculated values 
of R and the measured values of R still remains. 


F. In the search for commercially valuable minerals in Northern Ontario and 
Quebec, both single drill holes and drill holes in pairs have been used exten- 
sively for the purpose of placing the source internally in the medium or of 
placing both source and sink internally in the medium. Numerous measure- 
ments of the distribution of potential, of the interbow] resistance R, and of 
the resistivity p have been made experimentally, and in many cases the 
measured results have been compared with the values of R calculated from 
a suitably chosen initial value of p. 


1. The results of exploring the vertical drill hole No. U203 when C, is left 
at the bottom of the drill hole, P:, Pz moved up the drill hole, and the current 
is sent symmetrically in all directions about the drill hole are shown on Fig. 40. 


It is to be noted from Fig. 4b that in measuring the interbow] resistance, R, 
it is desirable to make the width of the interbow] cross section, viz., ) — a, as 
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small as possible, and the consequent interbowl differences in potentials are 
therefore also small in magnitude and the instrument which is used for mea- 
suring the interbowl potential differences must be correspondingly sensitive 
and accurate. 


D.D.H.*203 


eee ewes woees 


per cm. cube. 





Resistivities in 104 ohms 


100 200 300 400 500 600 700 


Distance from C, in feet. - C, down D.D.H*203 
835 ft. from collar ~ collar 1125 ft. below surface 


oe----0 Resistivities measured using potential electrode spacing of 50ft. 
Oe nn " " “ ” ” ” " 200 Ft, 


Fic. 4b. Resistivities from electrical logging. 


2. The layout for the investigation of the case where the source C; is placed 
in one drill hole and the sink C2 is placed in a neighboring drill hole is shown 
in Fig. 5. The vertical drill hole U203 shown in this figure is the same as 
that shown in Fig. 4a for which the ‘“‘logging’’ results are shown in Fig. 40. 
In the use of the pair of drill holes the source C; was placed in a fixed position 
near the bottom of drill hole No. U203, and the sink C2. was placed at succes- 
sively deeper positions in the adjacent drill hole No. U218. Potential contacts 
were placed in the drill holes above and below C; and C, in order to form 
equipotential surfaces at a short distance from C, and C; that become thus 
approximately spherical in shape. The potential difference and resistance, R, 
between these equipotential surfaces for successive positions of Cz were then 
measured. 


As has been presented previously, the experimental values of the interbowl 
resistance, R, at suitably chosen positions of C; and C2 were used to obtain the 
initial value of p which served as the foundation for the subsequently calculated 
values of R. An illustration of the results which were obtained for the layout 
which is shown in Fig. 5 is presented in Fig. 6, Curves I and II, for the 
experimentally measured values and also for the calculated values of R. The 
initial value of p was calculated from the measured value of R which was 
obtained, with the source C,; at the bottom of hole U203 and the sink at a 





A oC TAI 
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suitably chosen point in drill hole U218, using the modified formula (1), 
Section B. Subsequent values of R were calculated from this formula on the 
assumption of a homogeneous medium with the chosen initial value of p. The 
drill holes were not quite coplanar. Curve III shows the corrected values 
of R. 


To Potentiometer 






Genérator 
D.DH.*203 





P,, Py, P2,P2- Potential electrod 
C,,C, Current electrodes 


Fic. 5. Set-up with both source and sink in the medium. 


Even though the formula used to obtain the calculated values of R and p 
shown on Curve I was that developed for a homogeneous medium and therefore 
inadequate, it is quite apparent that a very informing comparison may be 
made of the experimentally measured values of R, as shown on Curve II, with 
the calculated values as shown in Curve I in which of course the experimentally 
measured values differed from the calculated values because of the variations 
from homogeneity in the medium. In particular, attention may be called to 


(i) The lowering of the interbowl resistance, R, with depth which from 
subsequent mining results was apparently due to disseminated pyrite. 


(ii) The pronounced evidence of the existence of dikes of high resistivity 
near the bottom of the vertical drill hole which also was found to be 
in accord with the results of drilling and in accord with the results of 
logging U203 as shown in Fig. 4d. 


(iii) Evidence of a test of the accuracy of the results of angling of drill holes 
is also shown in the position of the minimum value of the interbowl 
resistance, R, as presented in the Curve I of the values of this quantity 





. 


} 





j 
| 
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which is in almost precise agreement with the position of the experi- 
mentally determined value of the minimum resistance R, as presented 
in Curve II. The uncertainty introduced in the calculation of the 
resistance R by the assumption of a homogeneous medium was not 
very pronounced since the drill holes indicated in Fig. 5 commenced 
at the 1100 ft. level and were several hundred feet below the wet highly 
conducting overburden. In this case, then, the departures from homo- 
geneity as indicated by a comparison of Curves I and II, Fig. 6, were 
very closely in accord with the actual rock conditions indicated by the 
cores from these two drill holes. 








6000 1 ——- r 1 T 
| 

















t 
2 


Collars of drill holes 203 and 218 
were 1125 feet below the surface 
of the ground. Depth of wet 

overburden was approx. 100 fect 


Resistance in ohms from P,P; to P,P. 








eee 


{00 600 800 1000 1200 1400 
Distance in feet from collar of D.DH. U2I8 toCe 


Fic. 6. Inter-drill hole resistance between D.D.H. U-208 and D.D.H. U-218; source and sink 
in the medium. Curve I. The experimental values obtained using P, and P. as potential elec- 
trodes. Curve II. Theoretical values obtained by assuming the homogeneous medium of resistivity, 
p = 1.12 X 10° ohms per cm. cube. Curve III. Corrected values obtainable with the drill holes in 
the same plane. 








G. In many cases, however, drill holes which are used to provide the location 
for the source C,; and the sink Cz commence at the surface of the glacial over- 
burden. Now, where the layer of overburden is wet and highly conductive 
and covering a highly resisting layer of rock, it appears desirable to develop 
more precise formulas for the calculated resistivity p of a layered medium. 
The great difference between the conductivity of an overburden of water or of 
water-soaked glacial or alluvial material and that of a subsurface rock, even 
where the rock is wet, makes it especially desirable to make use of suitable 
electrical resistivity methods of exploration of the subsurface conditions. If 
the information which is obtained by the use of electrical layouts is correlated 
with information which may be obtained from other geophysical methods of 
exploration, e.g., magnetic, gravimetric, or seismic methods, together with 
information obtained from the cores, the location of the places where drilling 





= 
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is warranted may be made with confidence. Each drill hole should then be 
not only electrically ‘‘logged”’ but also electrically explored in every direction 
in the adjacent material with the source C; placed at various levels in the drill 
hole and the equal sinks placed on the surface symmetrically about C}. 


As the drilling proceeds, both source C; and sink C2 may be placed at various 
levels in each pair of drill holes and the calculated and experimentally measured 
values of R and p compared and thus extensive information may be obtained 
on the inter drill hole materials. 


H. \f the subsurface materials exist in an approximately layered structure it 
becomes desirable as an initial step to develop a suitable formula for use in the 
determination of R and p that will replace the formula that had been developed 
for use in the exploration of a homogeneous medium and which was used in 
obtaining Curve II, Fig. 6, viz., 


Pe agate (1-4 Se “oe I ). 
4n \a 6b L-a L-—-b 
The improved formulas which are suitable for use in the exploration of a 
medium of two horizontal layers of resistivities p; and ps, where the source is 
placed at various depths in a drill hole and the equal sink is placed on the 
surface or in a neighboring drill hole, were obtained by Rostocker. These 
formulas were subsequently simplified and made applicable in wider scope 
in the comparison of the calculated values of R and p with their experimentally 
measured values by Bernholtz for the case of two horizontal layers of resistivity 





. ° . ° ei 
p; and ps in which various values of the ratio — are assumed. 


P2 
The formulas which were obtained were then applied to the cases where 
a Po gens BN eee 
— has the values—— , —,—, and —. This provided a comparison of the 
ae 100’ 10° 5 2 


values of R and p of a layered medium; with the values of R and p of a homo- 
geneous medium. This comparison throws into relief the effect on the inter- 
bowl resistance, R, of the presence of a horizontal conducting layer or over- 
burden superposed on a lower layer of material of higher resistivity. This 
development is presented in Part II. 


Part II 


On the Distribution of Potential in a Medium of Parallel Layers 
with Source and Sink in Different Layers 


A. The problem considered here is that of obtaining more precise formulas for 
the potential distribution in a layered medium due to the presence of electrical 
sources and sinks placed in the medium. These formulas are used in calcu- 
lating the theoretical curves of interbowl resistance from which the inter- 
pretation of field results is made. In obtaining the potential functions it is 
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usual to assume the resistivity, p, to be a known function of position. The 
potential functions must of course satisfy the laws of electrodynamics and the 
prescribed boundary conditions. Here, horizontally stratified media will be 
considered in particular. 


B. The Method of Images 

Extensions of the method of electrical images originated by Maxwell have 
been found well suited to this type of problem. The method is found in the 
literature in various forms. Potential functions were obtained by E. Lancaster- 
Jones (8) in the case of the two layer earth for the Wenner electrode configur- 
ation. Explicit solutions for a three layer earth with a point source in the first 
layer were obtained by J. N. Hummel (3, 4) in the form of doubly infinite series. 
The method could at least in principle be extended to any number of layers 
of arbitrary and different thickness; however, beyond two layers the solutions 
become extremely complicated. D. O. Ehrenburg and R. J. Watson (1), 
making use of the image method in a different form, obtained solutions in the 
case of m layers where all layers are of equal thickness. 


The essential features of the method are readily revealed from a few simple 


cases. 





Fic. 7. Image for two layers (semi-infinite). 


1." In the case of two semi-infinite layers of finite resistivities and a source in 
one layer as in Fig. 7 the requirements of the potential functions Vo, Vi are 
(i) Vo, Vi must satisfy V?V = 
(ii) V; must have a simple pole at the origin. 
(iii) Vo must vanish at zg = + o-, 
V, must vanish at z = — ~, 


(iv) at s = hand for all values of r — ——- =— —_,, Vo= JV. 


It can be readily shown that solutions of the form 
pil 1—k s pol | k 

Vy=- <== = (a) 
4a Veit? + 22 2? 4n | V7? +2 Vr + (z- 2h)? 


NOTE: The lower case k in the text refers to the capital K in the figures. 








Vo= 
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1 av 
sp on 





satisfy the first three conditions and that J = — | dS over any sur- 


face enclosing the source. Condition (iv) requires that at z = h, 


Be cremate Re tee egies eee 5 


4a VP +)? 40 VP + h? pot pi 





Then — ——- = — — at z = h identically and (a) gives the unique solu- 
p2 OZ pi Oz 

tion of the problem. The important feature is that a singularity has appeared 
in V; at the point ry = 0, z = 2h. The boundary conditions could have been 
satisfied by assuming another source at this point of magnitude kJ. Then Vo 
could have been obtained from the original source with a reduction factor 
1—k. The problem is evidently analogous to the situation in geometric 
optics where the boundary plane is a surface of reflection factor k and trans- 
mission factor 1 — k. The electrical source and image are analogous to a light 
source and its image. However in this case the possible values of k are 
1 2 k 2 — 1, whereas in geometric optics k < 1. 


2. On investigation of the case of several layers as in Fig. 8, it is found that 
to satisfy the boundary conditions at plane (2) one must assume an image of 
magnitude ke/ as in Fig. 8a (8). But then the boundary conditions at plane (1) 
are no longer satisfied. To satisfy them one must assume an image source 
of magnitude k,kel as in Fig. 8(b) and then the conditions at plane (2) to be 
satisfied require another image k3kiJ and so on. In this way one can give 
significance to an infinite series of image sources and so,obtain the potential 
functions. 





Fic. 8. Development of images for three layers. 
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In general, then, to obtain the potential functions valid in a layer of resis- 
tivity p; it is only necessary to sum over all sources both real and image which 
an observer in the layer would see in the optical analogy with the appropriate 
transmission factors. 
pale zs J 

fj 


(8) 


i.e. Vz= 


where 
T;= transmission factor appropriate to the 7th layer, 
I;= magnitude of source, real or image, 
;= distance from a point in the ith layer to the jth source. 


C. Direct Solution of Layer Problems 

Alternately one may proceed in a straightforward manner to obtain solutions 
of Laplace’s equations which satisfy the boundary conditions even in the more 
complicated cases. This leads at once to cylindrical harmonic solutions due 
to the cylindrical symmetry associated with the special problem referred to in 
the introduction, and because of the nature of the boundary conditions leads 
also to Fourier integrals of Bessel functions. The resulting integrals may be 
evaluated by expressing them as a convergent series. This procedure gives 
the identical results obtained by the use of the image theory, as will be shown. 


1. In particular this straightforward method is employed here to obtain 
potential functions for a two-layered earth, each layer being homogeneous and 


Air P,=00 


Air-Earth Interface 







Pi,Vi 


Layer Interface 


P2,Ve 


Fic. 9. Two-layered earth—source or sink Fic. 10. Two-layered earth—source or sink 
in upper layer. in lower layer. 






isotropic and of constant resistivity (17). The source and sink may be in either 
layer (See Figs. 9 and 10). The prescribed conditions of this problem are, using 
cylindrical co-ordinates'7, g, z, as follows: 

ek 


r or 





(i) Vi, Vo, 





(ii) Vi>o as —&= Vaz 2 >O0, 
(iii) Va>O as E= VP +250, 
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: V; : : 
(iv) ——- = Oat z = + H, i.e., no outward normal current flow into the 


= 
~ 






air represented by a semi-infinite homogeneous layer of resistivity 





(vi) Vi= V2 at s = H —h, i.e., the potential functions are continuous. 





For (i) we assume a solution of the type V = e*“R(r) and obtain for R the 
equation R” ace: R’ + *R = 0, so that R = fJo(dr) | The solution Yo(Ar) 
’ | Volar) f 


= 0 and all values of sz. Hence V = e**J (Ar) and since X is unrestricted 


| 
must be rejected as it is infinite for r = 0 and this would make V infinite for 
because V vanishes only at ~ the most general solution is of the form 


v= |"Soaye™* + voret joa. 








» | 
Recalling the form of solution for a homogeneous infinite medium 
& 
» pt 1 
4m V/ 72 + 2? 
the general solutions appropriate to each layer are arranged to include this 
solution and perturbations. ; 
a i 
, 1 ‘ —s wi re 
Yyi=—|] — - + 1) Ai(A)e + Vilr)e f Jo(dr)dr ; (1) 
4 Vr + 2 0 { 


O<zsQt+H 


Vi= ao +{/ {a(rje* + yi(Aje—™} Joowyan | (2) 
T r- a 


H-h&2z§0. 





| 1 hos ‘ 
y= las ade +]. y0e*soOr)4 |, —o £2 H—h. (3) 


(1), (2) give the potential in the upper layer, and (3) gives the potential in the 
lower layer, when the source (or sink) is in the upper layer (Fig. 9). 





»T 1 se 
y= ; (de + wa(re*} eT ee 
1 ie | +e = ape I, {As(A)e + wW3(A)e | Joona git h RN AX Hy. (4) 


pol 


Ar 


4 1 i 
V2 = Iwas +]. Henan], — < z g AH —h. (5) 


(4) gives the potential in the upper layer, and (5) gives the potential in the 
lower layer, when the source (or sink) is in the lower layer (Fig. 10). 








- 
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Formulas (1) to (5) apply of course only to the case of the existence of either 
a source or sink only. 


The factor e~™ is excluded from V2 and V2’ because of condition (iii). These 


results are put in a more convenient form for applying the boundary condi- 
tions by making use of Weber’s integral. 


PP igs 1 
i. é 8 7,(\r)dy ” afrae » & 2 0 


r+? 


‘i " g 
and | e~*Jo(Ar)dd = -|- Ne Te(Ar)GA = eee: 
Zz J0 CTE +2 23 
The expressions for V1, V2, Vi’, V2’, now satisfy conditions (i), (ii), (iii), and 
application of (iv), (v), (vi) results in algebraic equations for 6:(A), ¥i(A), 
W2(A), O3(A), w3(A), Ya(A), and solving these equations we have: 
la. Source (or sink) in upper layer (Fig. 9). 
(i) Potential in upper layer: 
—— ell ae (ken. +o) sy, 
4In LVrt+e2 Jo 1 — ke? 


oO —2x —2r 
1 


(Ar)dy 


oe eJa0nr)dr], k= ps a 


l se he~2™ po + pi 
+H. (1’) 


ke~?*(1 +e") 
_— = =- + 5 catiniecai dl 
Vr+2 LC ee a 


) 
0<z 


M To(dr)dd 


j 
vim ST 
E 


—2dh 2\H 
+ < on e™Ioiwr)dr}; 
H-h€2€0. (2’) 


(ii) Potential in lower layer: 


y I 1 CO, —2\H —2dh : ae 
a 
y+ s 0 7 


—~a isi n—k (3’) 





1b. Source (or sink) in lower layer (Fig. 10). 


(i) Potential in upper layer: 


a of 1 —2dh 
Vy = 2 oT =+- ke) 7.09) dd 
0 











V7? +: + 2 j-- io 
(1—k)e™* ,, p2 — Pi 
- oe e*Jo(Ar)dX |, k = ———,,, 
0 — ; po t+ pi 
H-h2<dZ. (4’) 
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(ii) Potential in lower layer: 

ay pol 1 Me DH — ke ree) ‘ 

Y= i garst |, ger ee 
€z<LH-h. (5’) 





— © 


2a. It has been shown that (1’)—(5’) may be written as series by Stefanesco 
and Schlumberger. Then 




















Vi = 
ol 1 l = k a: e 
—} ————— + - + DS ——_ + 2 — 
4 lz +2 VP+(z— 2H)? *= 1/2 - + (2+2nh)? *=!1V r? + (2nh—2)? 
co kn oo Rk” 
$5 erences eeennngean eae Sepia |, 6 
n= 1/7? - + “(2nh -s Z- ~ 2H)? n=1v/ 72 + (2nh — 2 — 2H)? (6) 
F l 1 
gt fh Ua ea 
4 ( VP? + 2? VP + (2H — 2)? 
J = 1 1 
sea» E Kel Sciacca i. 7 
4a n= V r+ (2nh + 2H — 2)? VPr+ (2nh — 2)? (7) 





I oe 1 
v= - bo wi jeinia atlas cine | 
4 VP Toe + (2H — 2)? 








> oo bn pn 
pe Nias tceas 4 k : ; 
4r n=1 VP + (2H + 2nh — 2)? Vr + (Qnh 4 2)? (8) 
V's ” oe E a : & ee + ee 
| 44LVete2 Vr+ 2H -2h—2)? Vr+(2H — 2)? 
© La 
BE =| ‘ 
n=1V/ 72 + (2H + 2nh — 2)? (9) 


2b. These are the series obtained by J. N. Hummel (3, 4) by the image method. 
On examination of the above series for simple poles, image sources reveal them- 
selves on the r = 0 axis in the positions indicated in Figs. 11, 12. These are 
exactly the image sources predicted by the optical analogy and are of the 


required magnitudes. 
By performing the summations indicated in (8), Section B, 


f © n © n oO kn © n 
pr J RI kl I eT and 

















Wnetlis ET, + 
; 4a (nao PF" x Zs PPon-1 220 n= 1 PP’ on41 
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where PP's, =V r+ (Qnh — 2)? 
PP'ingy =VP + (2 + 2nh — 2H)? 
PP» =Vr? + (2 + 2nh)? 
PP. =V Pr + (2nh + 2H — 2)?. 





(0,0, 2H+2nh) 






(K"1)Pin+1 @ (0,0, 2nh+2H) 
(K"1)Pan (0,0, 2nh) 


(0,0, 2H+4h) 







(KI) Ps @(0,0,(2h+ 2H) 
(KI) P2 @ (0,0, 2h) 
(1) Py 0(0,0.2H y= 00 





(0,0, 2H+2h) 
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Fy 







(0,0.2H) _ 









(1) Qo ¢ (0,0.0) 


(KI)P, $ (0,0,-2h+2H) P,- 
(KI) Pe @ (0,0,-2h) 
woatae 





(0,0,-2h) 
(0,0,-4h) 






(0,0.-4h) (KI) Q2 
(K?1) Q, 


(K"1) P, n+i@ (0,0,-2nh+ 2H) 


(K®1)Pzn ¢ (0,0,-2nh) (KPI) Qon 


Fic. 11. Images for two-layered earth. Fic. 12. Images for two-layered earth. 


(K21) Py 












(0,0,-2 nh) 








These are the same expressions for V; and V2 previously obtained (6), (7) with 
the terms for » = 0 separated out. Further 


» _ pll—k) J SP kl kvl 
y= ey | t+ 
4r ] n=0 L PQaen PQon41 














vi EE Be 4 ae ee 
4x (PQ. PO: = POasss 
where PQ2, = r? + (¢ + 2nh)? PQon+1 = Vr + (2H + 2nh — 2)? 
PQ =V P+ 2 PO’, =VPr + (2H —2)?. 
(k has been preserved in the form k = ee as in the previous problem 
which accounts for the —ve sign in the 2nd term of V2’ and transmission factors 
(1 — k).) ° 


The series for V,/ and V2’ obtained in 2a, (8), (9) are evidently reproduced 
by the image theory. In (8), (9) the terms for » = 0 are separated out. 


D. Electrical Coring or Logging 
C. and M. Schlumberger and E. C. Leonardon (16) have applied the electrode 
configuration illustrated in Fig. 13 to study the formations traversed by a dr il 
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hole and have given the name electrical coring, later called electrical logging, 
to this process, which, in addition to being more rapid and economical than 
mechanical coring, can provide information concerning subsurface structure 
beyond that obtained from the examination of the drill core, since the electrical 
measurements are influenced by the surrounding media. 


In the case of impervious rocks where absorbed water is more or less sta- 
tionary the resistivity is a characteristic of the rock itself. However, in the 
case of porous formations such as sand, etc., where absorbed water circulates, 
resistivity is more a characteristic of the absorbed water, which will be high if 
the water is pure and low if the water has become an electrolyte from highly 


dissociated mineral solutes. 








Fic. 13. Schlumberger scheme for electrical coring. 
Ss 


In the case where the absorbed liquid is not a conductive electrolyte but 
instead an insulating medium such as, for example, porous formations with oil 
or natural gas, the resistivity is very high which is of some consequence in 
their detection by this method. 


In Fig. 13 if a, b >> the diameter of the drill hole the situation is ap- 
proximately equivalent to the case of a point electrode in the medium. On the 
assumption of a homogeneous semi-infinite earth, and neglecting the effects 
of the other electrodes at the surface and the air layer, this is a good approxi- 
mation. If measurements are made at considerable depth, we have radial 
flow, equipotential bowls as indicated in Fig. 13, and the applicable poten- 


tial function V = — —._ For the indicated measurements of potential dif- 
ference ; 
Va — Vi 4xab 


a ee ’ 
I h-—a 





(10) 


which is the Schlumberger formula for this set-up. 





| 
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As the three electrodes are moved up and down, holding a, 6 constant one 
would obtain always the same value for p computed from the observed values 
of the current and voltage if the medium were homogeneous (providing previous 
approximation conditions were fulfilled). However, when the medium is not 
homogeneous (this is the case of special interest) it is customary to make use 
of (10) to compute a quantity which is termed p,, the apparent average resis- 
tivity. Moreover, it is customary to attribute this characteristic quantity pa 
to only the material in the vicinity of and in particular between the two equi- 
potential bowls of radii a and 6 of the homogeneous earth. The question then 
is to what extent is this presentation justified? The change in p, will depend 
upon changes V, and V, from the case of a homogeneous earth (presuming 
the current J is held constant). It isa known result that V, and V, are affected 
by even remote inhomogeneities (such as a distant layer) to a considerable 
extent. In view of this the customary procedure of designating pg as the 
average resistivity of the material between the two spheres would appear 


unsatisfactory. 





Fic. 14. Two-layered case in electrical coring. 


However, it may be shown that the changes in V,, V, from the case of the 
homogeneous earth do actually occur in the vicinity of changes in p in such 
a manner that p,q will follow the changes in p, and significance may be given 
to pa as the average resistivity in the region between the equipotential bowls 
of radii a and b. 


With reference to Fig. 14 the potential functions are referred to the indi- 
cated co-ordinate axes and are readily obtained from the previous solutions 
by taking the proper combinations of formulas (6)-(9). The measurements 
of potential difference along the r = 0 axis permit the calculation of the 
apparent resistivity of the material adjacent to this axis. Potential curves 
have been computed for a particular case in Fig. 15 and comparison is made 
with the potential function for a homogeneous earth computed from 


08S) Se ae 
ale) 4r|2H—2| [2 
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where both electrodes and the image of the lower electrode are considered. 
Whereas there is a discontinuous change in the resistivity at the boundary the 
potential function remains continuous but changes in such a manner as to 
reveal the approximate position of the boundary. 


Potential 
-1.00 -.80 -.60 -.40 -.20 
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P, = 1000 meter A 


Depth in| Meters 





fo wee---- homogeneous 
; two- layered earth 
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f 2601. 
- 

a 





Les a, a 


Fic. 15. Potential distribution for a two-layered earth. 





E. Conclusions 

Explicit solutions for the potential distribution in a two layer earth have 
been obtained with sources or sinks in all possible positions. In the course of 
this development it has been demonstrated how the formal solution of Laplace’s 
equation leads to the results of the image method which reveals an isomorphism 
of the problem of electrical conduction in stratified media with a problem in 
geometric optics. 
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Part III 


Exploratory and Interpretative Procedure with Source and Sink 
Embedded in the Medium 


A. It is desirable to extend the electrical methods of investigation with the 
source and sink on the upper surface of a layered medium in modified form to 
the use of the drill holes by which the source and sink may be embedded in the 
medium. It is also necessary to develop the potential function in which the 
results of the experimental investigation may be incorporated and to derive 
a method of interpretation of the results by a comparison of the potential field 
measurement curves with suitable theoretical potential curves in order that 
the divergences may be correlated with the subsurface geological structure. 


B. An outline of the methods of investigation with the source and sink placed 
on the upper surface or embedded in the medium has been presented in Part I. 
The layouts are shown in Figs. 1, 2, and 3. 


A great deal of experimental and theoretical research has been carried out 
with the source and sink placed on the upper surface of a layered medium. 
The design of the necessary instruments ‘and the technique used in these 
researches are well established. Many variations in layout to meet field 


conditions have been devised. 


The extension of the subsurface use of sources and sinks involves no essential 
change in the design of instruments or technique. The essential objective in 
the interpretation of the results obtained is the determination of the subsurface 
geological materials and structure. For the subsurface exploratory work the 
problem of obtaining the curves of potential distribution and the methods of 
interpretation of the results follow a pattern similar to that for the case where 
the source and sink are placed on the surface. 


In both cases the formulas obtained for the potential function are complex 
even for the case of a geological structure as simple as that of two horizontal, 
homogeneous, and isotropic layers. In each case it is necessary to put the 
potential function in suitable form for numerical computation. Otherwise 
interpretation becomes very difficult. For surface exploration of a two-layered 
earth extensive numerical tables have been published by I. Roman (13). 
Graphical methods for systematic interpretation have been devised by G. F. 
Tagg (18, 19), with variations by I. Roman (12, 14), W. Longacre (9), and 
I. E. Rosenweig (15). 


For the extension of the resistivity methods of investigation to such cases 
as are shown in Fig. 5, i.e., with the source and sink embedded in the medium, 
no tables or collections of curves appear to be available for use in the inter- 
pretation of the results of measurements in the field. With a view to facili- 
tating such interpretation it appeared desirable to put the applicable potential 
functions in suitable form for calculation and curve tracing. 





ad 
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C. In this case, the assumed geological structure is a two-layered earth, each 
layer being homogeneous, isotropic, and of constant resistivity (see Figs. 9 and 
10). The distribution of potential in the medium is then given by formulas 
(1’)-(5’) of Part II. 


The formulas for the distribution of potential when the source is placed on 
the surface of the earth follows from formulas (1’)—(5’) by setting z = H = 0, 
with resulting simplifications in their form. The problem of deriving from 
these simplified formulas other expressions suitable for numerical calculations 
is correspondingly simplified, and as a result has been studied in great detail 
as stated above. Various rigorous expressions, either power series or asymptotic 
expansions, have been developed which are more or less suitable for calculations 
in different ranges of k. 


Application of the methods of simplification used in surface resistivity work, 
as given, for example, in Muskat’s paper (10), to formulas (1’)—(5’) of the 
subsurface case led to various expressions, not all, however, uniformly suitable 
for numerical treatment. For example, formulas (1’)—(5’) can be written as 
power series (see formulas (6)—(9).. These particular forms are not suit- 
able for numerical computations when k ~ 1, as the series then converge 
very slowly. It should be stressed, however, that for k small (~ 0) these 
formulas give a very convenient method for obtaining numerical results. Other 
expressions, suitable for numerical work when k ~ 1 can also be derived from 
formulas (1’)—(5’). For values of k between these, interpolated curves might 
give results sufficiently accurate for most field work. However, for the sys- 
tematic construction of extensive numerical tables the fact that the integrands 
in formulas (1’)—(5’) converge like negative exponentials (exp(— ax), a > 0) 
suggests that numerical integration would be the most suitable systematic 
method of obtaining results from these formulas. 


In any calculation the required formulas are selected from (1’)—(5’) and 
combined in the proper way. Thus the potential at any point is the algebraic 
sum of the potentials due to the source and sink. Since in practice it is poten- 
tial difference between two points that is measured, this is found by taking the 
algebraic difference of the potentials of the two points. 


D. A sample set of calculations was carried out first for the layout shown in 
Fig. 5, with C; and P, kept stationary in drill hole No. U203, Ci being 1950 ft. 
below the upper surface of the upper layer. C, and P: were moved down drill 
hole No. U218. The upper layer was assumed to be 100 ft. thick, i.e., h = 100 ft. 
The distance PiC;= P2C:,= 50 ft. This set-up corresponds to an actual case 
where the collars of drill holes U203 and U218 were 1100 ft. below the surface 
of the ground, and drill hole U203 was 850 ft. deep. These calculations were 
carried out for several values of k, in particular, 


(i) k = 0.9802 p3= 100 Pl, 
(it) 2 == OVSIR2 -. . p= 10 pi, 








GILCHRIST ET AL.: DISTRIBUTION OF POTENTIAL 23 


(ii) Rk = 0.6667... pa= 5, 
(iv) & = 0.3833: 2... 3 p=. S01, 
(v) R = 0.0000... p= pi. 


These calculations were repeated for the identical set-up as given above 
except that the drill holes were now assumed to begin at the surface of the 
upper layer. In addition, the calculations were carried out for pe= pi, using 
the approximate formula 


sisi aR ge PC Ply L= CiCs. b= CPi Cai. 
Pe ee 
a 6b L-a L-—b 
The results are plotted on Figs. 16 and 17 respectively. The graphs are for 
resistance against distance, since, if we divide both sides of any one of the 
formulas by the current 7, we obtain a resistance term on the left hand side. 
If the calculations are carried out for (1’)—(5’) in this form the significance of 
the results is independent of the values of the current J. In both cases for 
po p; a family of similar curves is obtained and the minimum points in all 
the curves occur for the same distance of C2 from the drill hole collar, namely, 
925 {t. This is to be expected since then P, is closest to Pz. The curves give 
some indication of the location of the interface between the two layers although 
it is possible to read off the depth to the interface approximately only. 


The curves of Fig. 16 show clearly the effect on the interbowl resistance, R, 
of the presence of a horizontal conducting layer or overburden superposed on 
a lower layer of material of higher resistivity. In particular, it may be noted 
that at depths well below the overburden the influence of the upper layer even 
when it is highly conducting is small. This implies that the simple formula 
for a homogeneous medium may be applied even when the medium is layered 
as assumed. This simplifies the computations involved, since the formula for 
the potential difference between source and sink in a homogeneous medium is 

5 2rR 
p= —$—$— —, which reduces to = Wager a when a = 


sa igh etee i Sapte 
a b L-a L-—-b a L-a 


P,C,= P2C. and L = C,C2, since b = L —a and L — b = a, from which 
numbers may be easily obtained. At depths approaching to, and on through 
the interface into the upper layer, the simple formula fails completely and the 
more precise formulas must be used. The necessary information on the values 
of the thickness and resistivity of the upper layer for the proper application 
of the complex formulas are obtainable from the cores and surface resistivity 
work. 


On Fig. 16 there is also plotted the actual curve obtained in the field from 
measurements which were taken at depths well below the surface of the earth. 
It is quite apparent that the same informing comparison may be made of the 
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experimentally measured values of R, with any of the calculated curves, in 


: ’ ; Ra (L— 
particular with the curve obtained from the simple formula p = sake \e 9) ’ 


as is done in Part I. L—2a 
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Fic. 16. Theoretical curve of interdrill hole resistance No. 1. 


E. As previously stated, the drill holes shown in Fig. 5 commenced at the 
1100 ft. level and were therefore several hundred feet below the wet, highly 
conducting overburden. In many cases, however, drill holes which are used 
to provide the location for the source C; and the sink Cz commence at the 
surface of the glacial overburden. This has been assumed to be the case in 
the second calculation the results of which are plotted on Fig. 17. Here again 
_ 2nRa (L—a) 


at depths well below the overburden the simple formula p rs 
— 2a 
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gives results in excellent agreement with the curves which represent the results 
obtained using the proper combinations of the exact formulas (1’)-(5’). 

The application of such theoretical calculations is sufficiently obvious not 
to require much comment. In the search for commercially valuable minerals 
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Fic. 17. Theoretical curve of interdrill hole resistance No. 2. 


the resistivity method has been widely used and its extension to subsurface 
work is now gaining ground. The success of this method depends on the 
existence of extensive tables and sets of graphs for purposes of interpretation. 
It is proposed to prepare such tables for at least certain wide classes of situ- 
ations which may arise in resistivity prospecting, and to apply these tables 
to show their usefulness in the interpretation of a considerable body of experi- 
mental interdrill hole results. 
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F, It has been presented that in the experimental electrical layout which is 
used to measure the interbow] resistance between bounding equipotential faces 
or bowls it is necessary to make provision for the following conditions: 


1. For a layout of the Wenner form or a symmetrically similar form the source 
and sink which are applied at the bounding equipotential faces should be equal 
in magnitude and preferably opposite in sign. 


2. For any layout the equipotential face or bowl of zero value should be an 
infinite plane face bisecting the region under exploration. If the ‘electrical 
current is supplied by a generator of low internal resistance and connected 
directly to the points of entrance of the current at C; and of emergence at C2 
and if the ‘‘stake resistance’ to earth of C; is approximately equal to the 
“stake resistance’’ at C; then the infinite plane equipotential face of zero value 
will bisect the region between C; and C;. The magnitude of the source C, 
and of the sink C; will also be equal. 


The necessity for making experimental provision for conditions 1 and 2 
becomes imperative, since the basis on which the simple formula for a homo- 
geneous medium and the more complex formula for a two-layered medium 
is the assumption that the source and sink are equal in magnitude and that 
symmetry of potential exists. 


For electrical exploration on the surface of the ground these conditions, 1 
and 2, may be provided for experimentally if a third point X, of entrance of 
current to the earth of low “‘stake resistance’’, be provided at some distance 
from the points C; and C, and placed symmetrically with reference to these 
points. Then the total resistance of the electrical circuit X + C, is made 
equal to the total resistance of the electrical circuit X + C2 by adjusting the 
“stake resistance’ at C; or at Co. Thus the “‘stake resistance’ at C, becomes 
equal to the “‘stake resistance’ at C2, and a state of electrical uniformity and 
symmetry is provided for a homogeneous medium. The experimentally mea- 
sured interbowl resistance between the bounding equipotential faces may then 
be compared with the calculated interbow] resistance for a single homogeneous 
or a two-layered medium with homogeneous layers which is obtained by the 


use of the formulas which have been developed. 


G. If drill holes are used to provide entrance and emergence of electrical 
current at subsurface points C; and C2, provision should also be made for the 
existence of the conditions given in F.1 and F.2. This may be provided by 
methods similar to those used when C; and C, were points on the uppermost 
surface of the ground. 


However, if the current and potential contact electrodes are of the form, 
size, symmetry, and distance of separation which are presented in Fig. 18 and 
which are shown in place in Fig. 5, the required conditions given in F.1 and 
F.2 will be provided for approximately, especially if the drill holes are filled 


with water. The experimentally measured interbowl resistance between the 
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bounding equipotential faces may then be compared with the interbow] resis- 
tance for a homogeneous or two-layered medium which is obtained by the use 
of the formulas which have been developed. If the medium is two-layered 


1. W- Weight to pull 
apparatus down pl 







| = 
drill hole ] j 
2. There are 3 brass strips | 
on each of the copper [L| IL 
tubes P,,C,,P,. S S 
| 30 feet 
| 
| | 
ode 
lit S S 
Lugs 
Copper tube cape 
3 feet 
Brass Strip 
mn Fa 
Ss T'S 


Fic. 18. Set-up for current and potential electrodes in drill holes. 


the simple formula is applicable only if the source and the sink are remote from 
the plane of separation of the layers. For positions of the source and the sink 
in the neighborhood of this plane, the more complex formula must be used. 
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ON THE ABSENCE OF PHOTONS AMONG THE DECAY PRODUCTS 
OF THE 2.2 MICROSECOND MESON: © 







By E. P. Hincxs ANpD B. PONTECORVO? 








Abstract 


An experiment is described which tests the hypothesis that the cosmic ray 
meson with a mean life of 2.2 usec. decays into an electron anda photon. Geiger 
counter trays are used to select mesons incident on a graphite block, and to detect 
decay products emerging from the graphite. The electronic circuits record 
delayed coincidences that correspond to a decay event occurring between 0.6 
and 5.3 wsec. after a meson is stopped. The absence of delayed coincidences of 
a type that could be attributed to the simultaneous emission of an electron and 
a photon, each of ~ 50 Mev., shows that the above hypothesis of the meson 
decay process is incorrect. The experiment also demonstrates the absence of a 
hypothetical unstable neutral meson among the decay products. 















Introduction 






The decay process of the ‘‘ordinary cosmic ray meson’’* has been the object 
of numerous investigations. The experimental study of the meson decay was 
encouraged by the theoretical work of Yukawa, and for many years seemed 
to give a brilliant confirmation of his views. According to Yukawa’s meson 
theory of nuclear forces, a charged particle of about 200 electron masses and 
integral spin is responsible for the intense short range forces that hold to- 
gether nucleons. When positive and negative mesons with masses close to 200 
m, were later found to constitute the hard component of the cosmic radiation, 
they were naturally assumed to be the predicted Yukawa particles. In ad- 
dition, according to Yukawa, the ordinary nuclear beta decay could be ex- 
plained by a typically quantum two-step process—virtual creation of meson 
by nucleons and subsequent radioactive decay of the ‘‘virtual meson’’—pro- 
vided the free meson is itself unstable against emission of an electron and a 
neutrino. 




















It was recognized that the mean life of a Yukawa particle that is required 
to explain quantitatively the beta process by the above mechanism is con- 
siderably smaller than the observed value—2.2 usec. (5, 23)—for cosmic ray 
mesons. Nevertheless, many reasons, among which were the very instability 
(29, 35) of the meson, its mass, and the fact that the charged particles resulting 
from the decay of the meson at rest have a penetration through matter that 









1 Manuscript received August 4, 1949. 
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Energy of the National Research Council of Canada. Issued as N.R.C. No. 2057. 
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Energy Research Establishment, Harwell, England. 

*There is strong evidence that this meson, which has a mass equal to— 200 electron masses and 
constitutes the bulk of the hard component of the cosmic radiation at sea level, is to be identified with 
the ‘‘u-meson”’ observed by Lattes, Occhialini, and Powell (Nature, 160: 453. 1947) among the 
products of disintegration of a heavier meson. In this paper the word meson will refer to the 200 
electron mass particle, unless otherwise stated. ' 
tSee, for example, reference (11). 
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cannot be very different from that of 50 Mev. electrons (6), supported the view 
that the meson decays into an electron and a neutrino, each of about 50 Mev., 
and is to be identified with the Yukawa particle. This picture of the meson 
decay process was generally accepted until, in 1947, Conversi, Pancini, and 
Piccioni (4) showed experimentally that the meson has an extremely weak 
interaction with nucleons and cannot be responsible for nuclear forces in the 
Yukawa sense. The strength of the interaction was found (8, 33, 34) to be 
~ 10-" times that required by the Yukawa theory. A consequence of the 
experiment of the Rome group is that the Yukawa interpretation of beta decay 
fails as well, since it would now predict much, too long beta lifetimes. The 
meson is then no longer the particle responsible for the nuclear beta process.* 
Consequently there is no reason to assume that charged mesons have integral 
spin (as the Yukawa explanation required), and, in particular, there is no strong 
reason for the assumption that the meson decays with emission of a beta par- 
ticle and a neutrino. 


The hypothesis may therefore be considered that the meson has a spin equal 
to 3(h/2mr). Sucha spin value would be required, for example, if the process of 
nuclear absorption of a single meson is accompanied by the emission of one 
neutrino in analogy with the process of electron capture by nuclei, as was first 
pointed out by one of the authors (28). This analogy between the processes 
of meson absorption and ordinary K-capture is suggested by the following 
facts. (i) The probability of capture of a bound negative meson (about 10° 
sec.', according to the experiment of Conversi and collaborators) is of the 


order of the probability of the ordinary K-capture process, when allowance 
is made for the difference in the disintegration energies and the difference in 
the volumes of the electron K-shell and of the meson orbit. (ii) Among a 
few pictures of a meson stopping in the gas of a cloud chamber, no star has 
been observed at the end of the meson track.f The absence of stars can be 
explained by the assumption (28) that a neutral particle, the neutrino, carries 
away most of the energy. 


A meson of spin 3(h/27) could, a priori, decay in a number of ways. One 
possibility would be the decay into one electron and two neutrinos.{ Another 
possibility would be that the 2.2 usec. decay process consists of the emission of 
one photon and one electron, each of about 50 Mev. according to the laws of 
conservation of energy and momentum. Radiative decay of mesons has been 
discussed by several authors (7, 9) from a theoretical point of view. 


It must be acknowledged that there is already some experimental evidence 


*The nuclear beta process can alternatively be described according to the original Fermi picture 
(without mesons). However, it must be noted that the failure of the Yukawa theory in explaining 
the beta process, which is being discussed here, concerns the particular meson of — 200 electron 
masses. It may well be that the Yukawa picture is correct, but that another type of meson is 
involved. 

tSee, for critical surveys, references (20, 25). 

TAt the time of writing this paper, the electron plus two neutrinos hypothesis, mentioned in 
1941 by Nordheim (24), is in favor, being strongly substantiated by the experimental observation 
of the form of the (continuous) energy spectrum of the decay electrons. 
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against the electron—photon hypothesis. The available analysis (3, 24) of the 
soft component of the cosmic radiation in equilibrium with its primary meson 
component suggests strongly that only a fraction of the meson energy is spent 
in ‘ionizing’ radiation (including y-radiation). However this argument, though 
strong, is indirect and, considering the uncertainties in our knowledge of meson 
physics and the interest of the problem, it was thought worthwhile to test in 
a straightforward way the electron—photon hypothesis. The experiment that 
was undertaken is reported in this paper; a preliminary note has already been 
published (13). Another experiment designed to detect photons in the meson 
decay has been published independently by Sard and Althaus (31, 32). 


In this paper it is assumed that the charged particle emitted in the meson 
decay is an electron. Recent experimental evidence obtained by the authors 
(16, 17, 18) after the completion of the work described here has confirmed the 
accuracy of this assumption. 


The Principle of the Experiment 


Two ways of designing an experiment capable of detecting photons in the 
meson decay have been considered. First, one may look for a delayed event 
that could only be due to the emission of the hypothetical photon after the 
stopping of a meson, without detecting the decay electron. Alternatively, one 
may look for the delayed coincidence due to the simultaneous emission of 
photon and electron in opposite directions. It is interesting to notice that 
the information given by the two methods is of somewhat different character, 
as will be explained later, and the two methods are complementary. 


The first method was selected by Sard and Althaus (31,32); the second 
method was chosen by the writers. From an instrumental point of view, it 
may be anticipated here that the first method requires an efficient anticoin- 


cidence system (implying a large number of counters), while the second method 


probably requires more electronic apparatus. ‘This consideration has influ- 
enced our choice. A more fundamental difference is discussed in the conclusion 
of this paper. 


After the method of electron—photon coincidences was selected, consideration 
was given to the most advantageous geometry. The essential condition to be 
achieved is that a vertical meson beam selected by a counter telescope be 
allowed to fall on an absorber whose minimum dimension is not large compared 
to the average range of the decay electrons. It should then be possible to detect 
electrons, and photons if they are present, emerging from opposite faces of the 
absorber. The two possibilities are illustrated in Fig. 1 (a) and (0). In (@) 
the meson is brought to rest in a ‘‘thick but narrow” absorber, while in (5) the 
meson is stopped in a “‘thin but wide” absorber. It turns out that arrangement 
(b) has a considerably greater ‘‘luminosity”’ than (a), essentially because of the 
divergence of the incoming meson beam. The arrangement illustrated in (0) 
has some unusual characteristics that will be discussed later. 
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ELECTRON 


(a) 


Fic. 1. Alternative orientations of meson absorber and decay particle counters for a vertically 
incident meson beam. The ‘‘thin but wide’ absorber arrangement (b) was selected for the present 
ex periment. 


The Experimental Arrangement 


The counter arrangement, which was placed in a concrete house in the Chalk 
River Laboratory (near sea level), is illustrated in the left part of Fig. 2. A 


Fic. 2. Arrangement of counters with simplified block diagram showing the essential functions 
of the electronic circuits. 


photograph of the apparatus is presented in Fig. 3. A 15 cm. thickness of lead 
above the apparatus has the function of absorbing the soft component and also 
of increasing (30) the number of slow mesons impinging on the meson telescope. 
The telescope is formed by the counter trays A and B. The mesons (about 
* 0.4% of the incoming beam) that are stopped in a graphite block 38 cm. by 
19 cm. by 8.5 gm. per cm.” thick, undergo decay. Graphite was selected 
because both positive and negative mesons decay when brought to rest in 
carbon (4). Decay electrons may be detected in either tray B or tray C. 
Decay photons, if present, will also be detected in either of these trays (B or C) 
whose efficiency for gama radiation is increased by introducing 2.1 mm. of lead 
between the graphite and both B and C. 
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A coincidence between A and B-—-which is referred to as an event (AB)— 
defines the passage of a meson. Coincidences between B and C delayed with 
respect to the event (AB) would result from the hypothetical electron—photon 
emission in the meson decay. The circuits were actually arranged so as to 
record the following delayed events. 


Fic. 3. Photograph showing the experimental arrangement with the counter trays A, B, and C 
in position. 


(B)ae1, or discharges of tray B occurring between 0.6 and 5.3 usec. after 
(AB). 
(C)ae1, or discharges of tray C occurring between 0.6 and 5.3 usec. after 
(AB). 
(BC)a.1, or coincidences of B and C occurring between 0.6 and 5.3 usec. 


after (AB). 


The function of the circuits is shown diagrammatically in relation to the 
counter arrangement in Fig. 2. It must be noticed that tray B has a twofold 
function: first, detection of the passage of a meson by a coincidence with A 
(event (AB)) and second, detection of a decay electron (or photon) following 
(AB). The use of tray B in this way was required by the adoption of the 
method (b) discussed in the previous section, giving high luminosity. Further- 
more, such an arrangement results in a considerable economy in the number 
of counters. A slight disadvantage is, of course, that one of the eight counters 
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of tray B, that through which the meson has passed, is insensitive to a decay 
particle because of the long dead time characteristic of Geiger counters. 


Electronic Circuits 


The electronic recording of (B)g.; events presents an interesting problem 
associated with the twofold function of tray B. In order to detect an event 
(B) 4.1 following a first event (B) by, let us say, 1 usec., it is clear that one could 
generate narrow (< 1 usec.) pulses and use circuits with a dead time smaller 
than 1 ywsec. While such circuits can be designed, a solution of the problem 
by the use of short dead-time circuits would have prevented the use of much 
standard apparatus available in the laboratory. A different and simpler solu- 
tion, which allows the use of long dead-time standard circuits, was adopted. 
The initial and the delayed pulse at the output of the Geiger-Miiller tray are 
distinguished as follows. In order that the first pulse should not trigger the 
input of the delay channel circuit, it is arranged that each delayed pulse is 
superposed upon the (properly shaped) initial pulse, and so rises to an ampli- 
tude double that of a single pulse. Amplitude discrimination between initial 
and delayed pulses becomes possible in this way. Each Geiger counter of 
tray B is connected to a single pentode, which is rapidly cut off by the begin- 
ning of a Geiger discharge. The square pulses of about 10 usec. duration which 
are thus produced at the plates are then mixed in a proportional stage. While 
the single pulse output, or ‘‘normal output’, is fed to the coincidence mixer 


detecting the event (AB), the ‘‘threshold output’’, at which the first pulse is 
biased off, feeds the delay channel. A full description of the circuit is presented 
elsewhere (18). 


A complete block diagram of the circuit functions is shown in Fig. 4. The 
amplifiers have a gain of about 1000, and saturate on all Geiger pulses. Their 
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Fic. 4. Detailed block diagram of circuits used to record the delayed coincidences. 
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purpose is to steepen the rise of the pulse so that the mixer inputs are triggered 
with negligible delay (({{ 0.1 usec.). 


The coincidence mixers are of fairly standard design, and are very similar 
to that described by Jelley (19). 


In the delay mixer, the circuit of which is given elsewhere (18), the input 
pulse corresponding to the event (AB) initiates a gating pulse 4.7 ysec. in 
width and delayed 0.6 usec. This gating pulse is then mixed separately with 
pulses from the B threshold output and from the C output, so that if a decay 
particle passes through B or C, between 0.6 and 5.3 usec. after an event (AB), 
a delayed coincidence (B),,; or (C) 4.) is registered. The delay mixer has, in fact, 
a third channel which was designed to register events (A)qe; from a threshold 
output in the A tray. Registration of the (A)4q.; events would have reduced 
the chance delayed rate by cancellation of all events (BC),4.; accompanied by 
(A)ae1, these (ABC),4.; events being produced by mesons discharging the three 
trays between 0.6 ysec. and 5.3 usec. after (AB).* However, since significant 
results were obtained without taking this step, no threshold output was added 
to A and the third channel was not used in the present experiment. 


The numerators consist of a mechanical register driven by a thyratron. 
Every delayed event is recorded in two separate numerators, one of which 
(only one is shown in Fig. 4,) supplies pulses to operate a pen recorder. 


The resolving times of the various mixers, and the initial delay and the 
gate width in the delay mixer, were measured with a calibrated oscilloscope, 
using a pulse generator which produced pairs of pulses whose time separation 
could be varied. The oscilloscope time calibration was checked by measuring 


TABLE I 


MEASURED VALUES OF THE MOST SIGNIFICANT TIME INTERVALS 


Position 1 | 0.6 psec. 


Initial delay of the J 
delay mixer \ Position 2 | 1.6 « 
Gate width of delay mixer | 1.7 5 
Resolving time of (AB) 0.7 7 
coincidence mixer 
Resolving time of (BC)del 7.0 2 


coincidence mixer 
the resolving time of one of the mixers, using random pulses from two counters 
placed far apart and counting 6-particles from two radioactive sources. In 
Table I are given the measured values of various significant time intervals 
characteristic of the apparatus. 


* The 1.5 cm. of lead between trays A and B was introduced to prevent a decay electron from dis- 
charging both B and A, at the time when it was planned to record (A )dei events, in addition to (B) det 
and (C)del. . 
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The Geiger Counters 


The Geiger-Miiller counters used in trays B and C (eight per tray) are all 
metal (copper cathodes) of the self-quenching type, 15 in. long (effective) and 


1 in. diameter. In tray A eight glass counters, self-quenching, 8 in. long and 
2 in. diameter were used because of their availability when the experiment was 
initiated. With their large diameter, the glass counters would not have been 
suitable for counting delayed events, on account of the counter time lags. For 


the detection of mesons, however, they were satisfactory. 


All counters are mounted inside aluminum boxes for mechanical support and 
electrical shielding. Fig. 5 is a photograph of one of the trays of copper 


counters with the lid removed. 


Fic. 5. View of eight copper counters mounted in aluminum box. (Tray C). 


The H.T. voltage, which is obtained from batteries, is applied to each counter 
wire through a separate resistor (107 ohms). In trays A and C the pulse from 
each counter is fed through a separate high voltage condenser (50 yut.) to the 
common output. In tray B, instead, each counter is connected by a condenser 
to the corresponding pulse-shaping tube of the ‘‘threshold amplifier’. 


Each counter was selected by measuring the plateau and by inspecting the 
pulse shape with an oscilloscope. For trays B and C a satisfactory plateau for 
the complete tray was achieved with a single voltage for the tray, but tray A 
was less satisfactory, and three different voltages were required. 


Counter lags, i.e., lags of the output pulses with respect to the passage of 
the ionizing particles, were investigated with a triggered oscilloscope. It was 
found that lags > 0.2 usec. were sufficiently rare to warrant the use of the 
counters for the present work. It should be noticed, moreover, that the most 
significant of the rates, the (BC)q.; rate, can hardly be influenced by lags even 
with a delay as short as that which was used. 
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Control Experiments 

The following controls were periodically made. 

The total counting rate for each tray, and also the rate of ‘‘triggering events’’, 
(AB), were measured twice a day without interrupting measurements. Typical 
values are given in Table II to give an idea of the intensities. 

TABLE II 


FREQUENTLY MEASURED PROMPT RATES (COUNTS PER MIN.) 


Meson absorber | (A) | (B) 


(C) | (AB) 


Graphite + lead between B and C | 1290 | 1000 
| 
| 


1230 360 


No graphite + lead between B and C | 1290 | 1215 | 1335 .}| 360 


Once a day the voltages on the counters were measured and the pulses from 
each counter tray inspected with an oscilloscope. 

On two occasions runs were intentionally taken with no voltage on the 
counters: no pulses were registered in -24 hr. 

The transmission of each branch of the circuits was checked occasionally in 
order to verify that the counting rates at the input and the output of the branch 
were the same. 

’ The significant resolving times, and especially the initial delay and the gate 
width in the delay mixer, were checked periodically. 

One experiment with an initial delay of 1.6 usec. (position 2 in the delay 
mixer), and with the usual gate width of 4.7 wsec., was performed as a check. 
The ratio of the net measured rate (B)ge. + (C)ae1, due to the graphite plus 
lead absorber, to that with 0.6 usec. delay, was consistent with a meson mean 
life of 2.2 usec. 

Sufficient controls were also available to avoid any possible misinterpre- 
tation of results due to power failures. 

The behavior of the threshold amplifier was checked periodically by mea- 
suring the bias curve (number versus bias) for the threshold output. This 
curve is discussed elsewhere (18) in detail. 

The constancy within the statistical error of rates (B)g.; and (C)q), as well 
as the distribution in time of the delayed events as shown by the records of the 
pen recorders, were good indications of the proper behavior of the apparatus. 


The Casual Rates 


The (B)q.; and (C)ge, rates not due to decay—the casual rates—are produced 
first by threshold amplifier lags* in tray B, and random counter lags, mainly 


* The term “‘lag”’ is used to refer to an accidental delay in the detection of a pulse due toa prompt 
event. Counter lags in tray B contribute little to the casual rate because at least one pulse must 
appear from B before a triggering event is established, and only a meson discharging two counters 
of B, one discharge occurring with a lag, could give a delayed count. ‘‘Threshold amplifier lags”’ 
may occur in cases in which the second B pulse has an abnormally low rising time. 
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in tray C, and second, by random particles discharging a tray between 0.6 and 
5.3 wsec. after the triggering event. The magnitude of the latter effect (the 
chance rate) can be calculated by a knowledge of the circuit characteristics, 
and of the various relevant undelayed rates, while the contribution from the 
first effect cannot be estimated a priori. 


However, the experiment when done without graphite plus lead between B 
and C should give the total (B)4.; and (C)qe; casual rates, except for a small 
contribution which arises from mesons stopped and decaying in the counter 
walls. By an analysis of all observations, and by computing the chance rates, 
the approximate values of the casual rates given in Table III were obtained, 
for the condition ‘‘without graphite plus lead’. The (B)ge; and (C)a, casual 
rates in the condition ‘‘with graphite plus lead”’ are expected to be 10 to 20% 
smaller than without graphite plus lead, because of the reduced prompt rates 
(B) and (C). The (BC),.; casual rate, however, is the same for the two con- 
ditions. 

TABLE IIT 
CASUAL RATES (COUNTS PER HOUR) 


Origin of casual counts (B)del (C)del | (BC)del 
| _ — 


Threshold amplifier lags 3.6 | 


Chance 


Counter lags | 1.0 3.1 
| 
| 


1.7 0.9 | 0.20 + 0.02 


Total casual rate D 4.0 | 0.20 + 0.02 
| 


The chance (B)q.; and (C)4.) rates are obtained from the following formulae: 
(Bygien* = [(AB) — (ABB)] [f (B)] r counts per hr. 
(cygie"** = [(AB) — (ABC)] (C) + counts per hr. 
where the prompt rates (Tables IJ and IV) are given in counts per min. and 7 
is the delay mixer gate width in seconds. The justification for the above for- 
mulae will be found, by analogy, from the considerations given below regarding 
the (BC){" rate. 


In Table III the casual (BC),.; rates have been included. These need some 
discussion as they are extremely important in the interpretation of the experi- 


TABLE IV 





PROMPT RATES (COUNTS PER MIN.) 


| 
(ABC) | 214 
(BC) 373 
| 


| 
(ABBC) | 14.5 


(ABB) 33.4 
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ment. First, it can be stated that there is no contribution to the casual 
(BC)aet rate due to counter lags and threshold amplifier lags, because the 
probability of simultaneous lags in C and B is negligible. As for the chance 
(BC) qe, it must be noted that (AB) events in which tray C is discharged simul- 
taneously (chiefly due to mesons passing through the three trays) can never 
be associated with a chance (C),4.; event, nor therefore with a chance (BC) qe 
event. This is so because of the long dead time associated with the C branch 
of the mixing circuit.* Then, to a first approximation, the number of chance 
(BC)a events is |(AB)—(ABC)] (BC) 7 counts per hr., where the prompt 
rates (counts per min.) are given in Tables II and IV. However, corrections 
must be introduced to the above expression for the following reasons. 













1. One of the counters of tray B is insensitive after having been discharged 
in the triggering event: a factor 7/8 is therefore included to express the 
fact that for a subsequent particle only seven counters are available. 






2. An oblique particle discharging in addition to A two counters of tray B 
(event (ABB)) will give both a normal and a threshold output, making the 
B branch of the delay mixer insensitive to a subsequent pulse. 







Taking into account these considerations, the chance (BC),4.; rate may be 






written 


(BC) iii" = [(AB) — (ABC) — (ABB) + (ABBC)] [3(BC)] + counts per hr. 






The relevant prompt rates not already given in Table II are given in Table IV; 

they do not change significantly with the insertion of the graphite plus lead 

between B and C. Evaluating the above expression we find finally: 
(BCG = (BC)R"~ = 0.20 counts per hr., 


del 








with an estimated error ~ + 10%. 








Results and Discussion 






Runs with and without the graphite plus lead absorber between B and C were 
alternated, and the measured counting rates with their standard deviations 






are given in Table V. 
The observed rate (BC),.; could be due to the following causes: 






(i) Genuine electron—photon coincidences from the meson decay, i.e., the 





effect looked for; 
(ii) Single decay electrons which traverse both trays B and C; 





(iii) Casual events. 





The casual (BC)4.) rate given in Table III is expected to be, to a very good 
approximation, independent of the presence or absence of graphite plus lead, 
as mentioned before. An inspection of Table V shows that while the observed 






* It is on account of this fact—the ‘‘anticoincidence function”’ of long dead-time circuits-—that 
an anticoincidence (AB — C) was not used as a triggering event. For the same reason the chance 
(B)del rate is greater than the chance (C)dei rate (Table ITI). 
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(BC) qe, rate with the graphite plus lead absorber is about equal to the estim- 
ated casual rate, the (BC)4.; rate without the absorber is substantially higher 
than the casual rate. This increase observed in the absence of absorber must 
be due to a delayed radiation considerably absorbed by the thickness of 
graphite plus lead used. To verify that single decay electrons (cause (ii)) are 
responsible for this effect a subsidiary experiment was performed. 


TABLE V 


DELAYED COINCIDENCE COUNTING RATES (COUNTS PER HOUR) 


= 
C)del (BC )del 


Meson absorber (B)del | (C)del | (B)det + ( 


With graphite plus lead | 
(104.2 hr. of observation)| 11.93 + 0.34 | 12.26+ 0.34 24.19+0.48 | 0.2140.05 


| 
Without graphite plus lead] | 
—(77.2 hr. of observation)| 6.48 + 0.29 64 + 0.2! | 


Net effect due to decay 
electrons from graphite] 5.45+0.45 | 7.62+0.- 13.07 + 0.62 
plus lead 


11.12 + 0.38 | 0.43 + 0.08 





Consider the counter arrangement (Fig. 2) without the graphite plus lead 
absorber. Decay electrons which discharge both B and C, contributing to 
(BC) ae, must be attributed to the decay of mesons which having traversed 
A and B are stopped in material below C, but do not discharge C.* Some 
decay electrons from such mesons will travel upwards through C and B. The 
number of mesons “‘leaking’’ through tray C without discharging it was in- 
creased by removing the voltage from some of the counters of tray C, while 
the fraction stopping just below C was increased by placing a thickness (8.5 
gm. per cm.*) of graphite below C. The delayed rates were measured and an 
analysis of the results shows that the excess (BC),4., rate noted in our experi- 
ment without graphite plus lead absorber may be attributed to this cause (ii). 
In the presence of graphite plus lead, the effect is not detected because, with 
the geometry used, the thickness of absorber between B and C is sufficient to 
reduce substantially (15) the number of decay electrons which can traverse 
both B and C. 


It is apparent, therefore, that the change in contribution of effect (ii) to the 
(BC) 4.) rate prevents an estimate of effect (i)—electron—photon coincidences 
being made by comparing the rates (BC) .; with and without graphite plus 
lead. Consequently it is necessary to compare the measured rate (BC) .; with 
graphite plus lead with the estimated casual rate, (BC){7{“". The differ- 
ence observed (BC): rate with graphite plus lead minus casual rate 

casual +. 
rie | + 0.06 


counts per hr. 
|— 0.01 


A = 0.01 


* Because of the long dead-time anticoincidence function of the C circuit. 
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A is to be interpreted as the observed contribution of the hypothetical decay 
into an electron and a photon, and must be compared with the expected effect 
from such a decay process. 


This effect can be estimated by taking as a starting point the observed 
numbers of decay particles (B)i{{, (C)i¢j, from graphite plus lead, which are 
given in the last line of Table V. Let & be the efficiency of one of the lead 
covered trays for gamma radiation of 50 Mev., and, assuming the photon hypo- 
thesis, let 8 be the probability that when a decay electron is detected in B the 
accompanying photon impinges on the lead covering C, and let y be the prob- 
ability that when a decay electron is detected in C the accompanying photon 
impinges on the lead covering B. The rate (BC),,.; due to electron—photon 
coincidences will be 

> = BE (B)5 + vé(C)5 counts per hr. 


Apart from a small correction due to the absorption of photons in graphite, 
6 and ¥ are essentially geometrical factors which were estimated to be approxi- 
mately 0.6 and 0.5 respectively. These figures represent averages over the 
absorber; the difference is due to the “‘dead counter”’ effect in B. 


The gamma efficiency — of a counter tray covered by 2.1 mm. lead for 50 
Mev. photons can be estimated in a straightforward way. Pair production is 
by far the main contribution to the total absorption process for such photons 
in, lead. Since the thicknesses of lead involved are small compared with the 
range of the most energetic electron (positive or negative) of the pair, the 
efficiency is only very slightly less than the probability that a photon produces 
a pair in the lead covering the tray. The absorption coefficient for a 50 Mev. 
photon in lead is taken to be 1 cm.~! (12), from which we find that £ is between 
15% and 20%. We can now arrive at the value > = 1 count per hr. with an 
estimated uncertainty of about + 30%. 


Conclusion 


A comparison of the ‘‘expected value”’ = of the rate (BC),,; with the observed 
value A leads to the conclusion that the 2.2 usec. meson decay process does not 
consist of the emission of an electron and a photon, each of about 50 Mev. according 
to the laws of conservation of energy and momentum. 


Our result refers to both positive and negative mesons, because in carbon 
mesons of both signs undergo decay. This is in agreement with the results of 
Sard and Althaus (31, 32) and of Piccioni (26, 27), which, however, refer only 
to positive mesons. 


The results of the present experiment and of the independent experiment 
of Sard and Althaus are also complementary. Substantially different gamma 
efficiencies and different “‘sources’’ of decay particles were used. In addition 
a rather fundamental difference should be emphasized. Sard and Althaus’ 
experiment is an attempt to detect delayed gamma radiation in the meson 
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decay process, without any correlation with the decay electron. In this sense 
their experiment is more general than the present experiment. On the other 
hand, the present experiment was designed with the intention of minimizing 
the chances of detection of bremsstrahlung from the decay electron ;* in fact 
bremsstrahlung photons would be emitted primarily in a forward direction 
by the electron while our experimental arrangement is sensitive only to simul- 
taneous emission of an electron and a photon in approximately opposite direc- 
tions. The negative result obtained in this research has permitted the authors 
to design another experiment to detect the bremsstrahlung from the decay 
particles. This will be reported in a separate paper (18). 


Another conclusion which can be reached from the results of this experiment 
concerns the stability of a hypothetical neutral meson. Several authors (2, 22) 
have ventured the hypothesis that the decay of the meson might involve the 
emission of a neutral meson of mass about 70 Mev. This was mainly suggested 
by two cloud chamber pictures (1, 2) of decay electrons whose energy was 
about 25 Mev. Greisen (10) has argued that “numerous phenomena which 
have defied even qualitative explanation up to the present are made under- 
standable in the light of the neutral meson hypothesis’, provided the neutral 
meson decays immediately into two photons. The gamma instability of a 
neutral meson has been predicted on theoretical grounds (7, 9). The experi- 
mental results reported here lead to the conclusion (14) that either the neutral 
meson is not emitted in the 2.2. usec. meson decay process, or, if it is emitted, 
it does not decay into two photons with a mean life < 107° sec. A similar 


conclusion, referring only to positive meson decay, was reached by Sard and 
Althaus (32). 


The absence of photons in the 2.2. usec. decay process is in agreement with 
the evidence arising from an analysis (3, 24) of the genetic relationship between 
hard and soft components of the cosmic radiation in the lower atmosphere. 
Actually it would have been difficult to interpret the cosmic ray results on this 


topic if the total energy of mesons decaying in flight were emitted in the form 
of ionizing radiation. 


Our result definitely eliminates one of the decay processes that could have 
been postulated for a meson of spin 4 (h/27). However, a half-integral meson spin 
still presents a considerable attraction as was pointed out in the introduction 
to this paper, and in more detail in reference (28). Another process involving 
such a value for the meson spin, and currently in favor, is disintegration into 
an electron and two neutrinos. Recent experimental evidence supporting this 
hypothesis includes (a) the establishment of the fact that the charged decay 
particles have electronic mass (16, 17, 18), and (0) the form of the energy 
spectrum of the decay electrons (21). 

* In the experiment of Sard and Althaus the bremsstrahlung was discriminated against by having 
an anticoincidence counter tray on the path of the photon which would be discharged by an accom- 


panying electron. However this discrimination is much less severe than the one used in the experi- 
ment described here. 
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HIGH TEMPERATURE CALORIMETRY 
I. A NEW ADIABATIC CALORIMETER' 


By L. D. ARMSTRONG? 


Abstract 


In this paper is described a new calorimeter for the measurement of specific 
heats at high temperatures, by the adiabatic method. The advantage is that 
the specific heat at a definite temperature can be determined by a measurement 
taken over a small temperature interval, with a precision of 1% or better, 
throughout the range 400°C. to 800°C. This permits a study of specific heat 
anomalies in this range. 


Introduction 


Two general methods have been used for the measurement of specific heats 
at high temperatures. One is by means of the adiabatic calorimeter, which 
has been highly developed by Moser (4). In this method the calorimeter 
containing the sample is heated electrically, and heat exchange with the sur- 
roundings is prevented by a jacket which is kept at the same temperature as 
the calorimeter. Some investigators (6) have modified this method by at- 
tempting to correct for the difference in temperature between the calorimeter 
and the surroundings, but only with limited success. Second is the so-called 
“drop method’’, in which a sample is heated to a known high temperature, 
and then dropped into an insulated calorimeter. The heat given off by the 
sample may be measured by the rise in temperature of the calorimeter, or an 
ice calorimeter may be used. The advantage of the adiabatic calorimeter is 
that it is possible to measure specific heats directly over small temperature 
intervals, and so to obtain the details of the specific heat — temperature curve. 
However, up to the present time accurate work with this type of calorimeter 
has been confined to temperatures below 500°C., because of problems of thermal 
insulation at higher temperatures. 


The high temperature adiabatic calorimeter to be described in this paper 
was constructed primarily to measure the atomic heats of the transition group 
of metals, from chromium to nickel. These metals are interesting for a number 
of reasons, but reliable high temperature data have hitherto been available 
only for iron and nickel. As the atomic heats are likely to be abnormal, differ- 
ential values are essential for theoretical discussion, and the adiabatic method 
is preferable in this case. By the use of techniques which were developed for 
measurement of thermal conductivities at high temperatures (2), it has been 


possible to make reliable determinations of specific heats up to 800°C. 


1 Manuscript received August 29, 1949. 
Contribution from the Department of Physics, University of Toronto, Toronto, Ont. 
Now at RCA Laboratories, Princeton, New Jersey, U.S.A. 
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Design of the Calorimeter 






The adiabatic calorimeter which has been constructed is shown in Fig. 1. 
The sample, A, could be either a machined cone of solid metal fitting closely 
into the tapered copper container, B, or it could be in the form of powder or 
loose pieces packed firmly together. In the latter case the container was filled 
with an inert gas, usually helium. The sample and container were heated by 
the calorimeter heating element, C, consisting of No. 32 Chromel A wire wound 
on pipe-clay tubing. The heater was insulated with alundum cement, and 
was surrounded by a thin copper tube fitting into a hole through the center 
The calorimeter proper, B, was completely surrounded by a 











of the sample. 
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Fic. 1. High temperature adiabatic calorimeter. 





copper radiation jacket, D, which was kept at as nearly as possible the same 
temperature as B by means of three heater elements, EFG. The jacket was 
supported by a perforated steel tube, J, and was enclosed in a furnace, K, 
which could be evacuated or filled with an inert gas, as desired. A thermo- 
couple type vacuum gauge was used to indicate the pressure in the furnace 






tube. 


In assembling the calorimeter, extensive use was made of a butt-welding 
apparatus similar to that described by Armstrong and Dauphinee (1). This 
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unit proved to be invaluable, and it would have been impossible to dismantle 
and reassemble the calorimeter to make adjustments or to change samples in 
a reasonable time without it. 


Temperature Measurements 


The temperature of the sample was determined by means of a platinum 
to 90Pt-10Rh thermocouple, attached to the lid of the calorimeter. The wire 
for this thermocouple was obtained from Johnson, Matthey and Company, 
and was guaranteed to better than 1° up to 1000°C. The e.m.f. was measured 
with a Leeds and Northrup type K2 potentiometer, and the e.m.f.—temperature 
tables of Roeser and Wensel (5) were used to find the temperature. Additional 
Pt to Pt-Rh thermocouples were attached to the inside of the guard jacket, D, 


To 
: Differential 
Selector Switch Couples 


POTENTIOMETER 


Fic. 2. Thermocouple circuits. 


and to the jacket lid. These served both to check the calorimeter thermo- 
couple and to indicate the temperature distribution over the jacket. A com- 
mon cold junction for the three thermocouples was kept in an ice bath. 


In order that the temperature of the guard jacket could be adjusted quickly 
and accurately, differential thermocouples were attached to the top, side, and 
bottom of the inner container, with their opposite junctions at corresponding 
points on the inside of the guard jacket. These couples were made from a 
gold—palladium alloy which gives a high thermal e.m.f. against copper. The 
copper container and jacket acted as a common return lead. 


The measuring circuit for the thermocouples is shown in Fig. 2. The selector 
switch was connected so that it also acted as a reversing switch for the hot 
junctions of the three Pt and Pt-Rh couples. By means of a double reversing 
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switch, the common cold junction and the supply battery for the potentiometer 
could be reversed simultaneously with the hot junction. If there were any 
e.m.f.’s in the circuit, outside the thermocouples and battery, this triple re- 
versal caused a deflection of the galvanometer after it had been once balanced. 
As it is very difficult to avoid stray thermal e.m.f.’s, these were compensated 
by a small e.m.f. introduced in the portion, A, of the circuit. The compensating 
e.m.f. was adjusted until revérsing the circuit caused no change in the balance 
point of the potentiometer. The compensation was then equal and opposite 
to the stray thermals, and the potentiometer reading gave the correct e.m.f. 
for the thermocouple. 


The differential thermocouples were connected directly to the galvanometer, 
so that a zero deflection indicated a correct balance between the temperatures 
of the calorimeter and the jacket. Compensating e.m.f.’s similar to that in 
the potentiometer circuit were introduced into one lead of each differential 
couple, to correct for any small differences between them. 


All thermocouple and heater wires were led to the bottom of the supporting 
column through pipe clay tubes, and were taken out of the furnace through a 
high temperature vacuum seal of the type described by Armstrong and 
Dauphinee (2). 


Control and Measurement of Heat Input 


' The desirable rate of heat input to the calorimeter depends on the speed and 
accuracy with which temperature equality between the calorimeter and jacket 


can be adjusted and maintained while the temperature is rising. It was decided 
to use throughout a heating rate of about | centigrade degree per minute, and 
the calorimeter and jacket heaters were designed accordingly. The calorimeter 
heater was operated from storage batteries, and the voltage and current were 
measured with a calibrated Weston voltmeter and a Cambridge potentiometer, 
respectively. The jacket heaters were operated from the 110-v. d-c. mains. 
All the heaters were controlled by series rheostats. 


The heating time was measured with a laboratory stop clock, which was 
checked regularly against a chronograph operated by a standard clock. The 
stop clock was found to be reliable to within } second in 10 min., as long 
as it was kept well wound. 


The furnace heating element was made from Chromel A ribbon, wound on 
asbestos insulation over the iron furnace tube. It was designed to give 1500 w. 
at 100 v., and was controlled by a Variac transformer with a series rheostat. 
This was found adequate to maintain operating temperatures up to 800°C. 
The furnace temperature was measured with a Chromel P-Alumel thermo- 
couple, resting against the inner side of the furnace tube. 


Method of Measurement 


The adiabatic method requires accurate adjustment of the various controls, 
and reading of the instruments, in rapid succession during any heating interval. 
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Moreover, it was desirable that it should be possible for one man to operate 
the apparatus. Consequently the controls and indicating instruments were 
laid out on a carefully designed panel, to give maximum convenience in the 
sequence of operations. It was found to be feasible for one man to carry out 
all the necessary operations, but only after considerable practice, and with 
some sacrifice in accuracy. Fortunately, assistance in the taking of readings 
was obtainable during the latter stages of the work. 


In preliminary trials it was found possible to keep the calorimeter and jacket 
at the same temperature to within 0.1°C., with a heating rate of 1° per minute. 
However, using a sample in the form of small pieces, and with an atmosphere 
of helium in the system, it required about four minutes for the temperature of 
the calorimeter to come to equilibrium after a heating period. For this and 
other reasons it was difficult to adjust the temperatures of the furnace and 
calorimeter so that no drift occurred when the calorimeter heaters were turned 
off. Consequently the practice was adopted of adjusting the different tem- 
peratures until the drift was less than 0.02° per minute. The remaining drift 
was then measured before and after the heating interval, and a correction 
applied. 


In view of the manner in which the calorimeter behaved, two different 
routines of measurement were adopted, and used alternately. In the first 
method a heating current was applied to the calorimeter for an interval of 
about 10 min., and the rise of temperature was determined. In this case each 


measurement involved the following steps, of which Items 3 to 9 had to be 


carried out in rapid succession. 


1. The furnace was stabilized at the required temperature. 

2. The calorimeter and guard temperatures were adjusted until a drift of not 
more than 0.1° occurred in five minutes. 
The temperature of the calorimeter was read accurately, using the pro- 
cedure of reversing the thermocouple. 

4. The sample heater was turned on, and simultaneously the time clock was 
started. 
The temperatures of the guard jacket were maintained the same as that 
of the calorimeter by adjusting the jacket heaters. 
Concurrently the applied voltage on the calorimeter heater was main- 
tained constant, and the voltage and current read. 
The calorimeter heater was turned off at the end of the required interval. 

8. The guard temperatures were kept equal to that of the calorimeter until 
the latter became constant. 
The final calorimeter temperature was read accurately. 
The drift in calorimeter temperature was again observed, for at least five 
minutes. 
The net temperature change was calculated, using corrections from Items 
2 and 10 if necessary. 









































American Chemical Society. 






3. 


}. The calorimeter heater was turned off. 





the form of synthetic sapphire. 
and has no known high temperature transitions. 
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In the second method, heat was supplied to the calorimeter at a constant 
rate, and the rate of rise of temperature was determined. 
3 to 7 of the procedure were replaced by the following. 


In this case Items 


The calorimeter heater was turned on, the voltage maintained constant, 


and the voltage and the current were read. 


The temperatures of the guard jacket were maintained the same as that 


of the calorimeter. 


When the rate of increase became steady the time required for the cal- 
orimeter temperature to change by a specified amount was observed. 











Calibration of the Empty Calorimeter 










Two series of measurements were made to determine the heat capacity of 
the empty calorimeter, one series with the calorimeter empty, and one with 
it filled with a standard material of known specific heat, aluminum oxide, in 
This material, if colorless, is of high purity 
The 


specific heat has been 
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carefully measured by Ginnings and Corruccini (3), and it has been 
mended as a standard material by the Conference on Calorimetry 
The results obtained are shown in Fig. 
it will be seen that the two methods agree within experimental error. 


Heat capacity of empty calorimeter. 






recom- 
of the 
3, and 

The 


smoothed curve was used to give the heat capacity of the calorimeter for later 






measurements with other samples. 


In drawing it the measurements using 
aluminum oxide were given slightly more weight than those with the empty 
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calorimeter. The reasons for this were twofold. First, operating conditions 
with aluminum oxide were similar to those with other samples in the calori- 
meter. Secondly, these measurements were made in the later stages of the 
work, when practice in operating the apparatus had considerably improved 
the technique. 


An interesting feature of the heat capacity for the empty calorimeter is the 
irregularity between 150° and 200°C. The cause of this has not been estab- 
lished, but, in the absence of other reasonable possibilities it is tentatively 
attributed to an anomaly in the specific heat of the alundum cement, in which 
the heater element is embedded. 


Experimental Errors 


From a consideration of the precision with which the temperatures and 
heating current could be measured, and with which the heat exchange could 
be eliminated, it was estimated that the maximum error of an individual 
specific heat measurement would be less than 2%, and that values derived 
from a smoothed specific heat curve should be correct to within 1%. The 
accuracy actually attained is indicated by the scattering of the experimental 
points with respect to the smooth curves, as well as by the test with synthetic 
sapphire, and appears to be somewhat better than the above estimate. The 
average deviation of individual determinations from the smooth curves is 
about 0.5% below 400°C., and increases to not more than 1% at 800°C. 
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HIGH TEMPERATURE CALORIMETRY 


Il. ATOMIC HEATS OF CHROMIUM, MANGANESE, AND COBALT BETWEEN 
0° AND 800°C.! 


By L. D. ARMSTRONG? AND H. GRAYSON-SMITH?® 


Abstract 


The atomic heats of chromium, manganese, and cobalt have been measured 
up to 800°C. by the adiabatic method described in Part I of this series. The 
specific heat of chromium is regular from 0° to 800°C. The measurements with 
manganese show the sharp a-@ transition at 717°C., with a latent heat of 450 
cal. per mole. Typical supercooling occurs on the reverse transition. Cobalt 
shows an anomaly due to a gradual transition, which is known to be a change in 
crystal form. The new data obtained permit a comparative discussion of the 
atomic heats of the transition metals, chromium to nickel. All these have values 
Cy> 3R at high temperatures. After allowing for the effects of ferromagnetism, 
the excess specific heats of cobalt and nickel are accounted for by the conduction 
electrons. This is not true for chromium and manganese, for which metals there 
must be some additional source of internal energy. It is tentatively suggested 
that these two metals may have antiferromagnetic transitions at temperatures 


above 800°C. 
Introduction 
The adiabatic calorimeter described in Part I of this series of papers has been 
used to determine the specific heats of the metals chromium, manganese, and 
cobalt up to 800°C. The adiabatic method has proved especially advantageous 
in the cases of manganese and cobalt, since it has given information on tran- 


sitions which take place in these metals. The metals used were the purest 
obtainable, since considerable doubt as to the validity of the measurements 
exists unless the purity and state of the metal are known. 


Chromium 


The chromium used was in the form of electrolytic flakes, obtained from 
Johnson, Matthey and Co. The stated maximum amount of impurity was 
0.1%, chiefly chromous oxide and adsorbed hydrogen, with some calcium and 
sodium. The atomic heats are plotted against temperature in Fig. 1, and 
smoothed values of the specific heat are given in Table I, column 2. 


The atomic heat of chromium is quite regular over the range of temperatures 
concerned, the principal point of interest being the steady increase at high 
temperatures, to values considerably greater than the classical theoretical 
value C,= 3R. Two previous series of measurements are available for com- 
parison. Wiist, Meuthen, and Diirrer (17) have made measurements by a 
drop method with an ice calorimeter, on metal of unspecified purity. Umino 
(16) has measured chromium of about 95% purity, also by a drop method. 

1 Manuscript received August 29, 1949. 
Contribution from the Department of Physics, University of Toronto, Toronto, Ont. 


Now at RCA Laboratories, Princeton, N.J. 
Now Professor of Physics, University of Alberta, Edmonton, Alberta. 
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ra 


Mean values‘obtained by these investigators are included in Table I, columns 
3and4. The agreement with the measurements of Umino is good, considering 
the impurity of the metal which he used. The discrepancy with the measure- 
ments of Wiist et al. may be due to impurities in the metal used by these in- 
vestigators. 
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Fic. 1. Atomic heat of chromium. 


TABLE I 


THE SPECIFIC HEATS (IN CAL. PER GM. PER °C.) OF 
CHROMIUM, MANGANESE, AND COBALT 


Temperatuie, | Chromium ‘Biapmpnete Cobalt 
°C, poe aie: 
A&GS| Umino | Wiist A &GS| Wiist | A&GS| Umino 


Wiist 


0 | 0.1038 | 0.1023 | 0.1111 | 0.1204 | 0.0983 | 


0.0912 


200 | 0.1216 | 0.121 | 1157 | 0.1274 .1305 | 0.1132 | 0.114 0.1073 


ce] 
100 0.1146 | 0.118 | 0.1090 | 0.1200 | 0.1254 | 0.1056 | 0.110 | 0.0993 
| 


300 | 0.1259 | 0.123 | 0.1224 | 0.1336 | 0.1356 | 0.1200 | 0.118 | 0.1154 

100 1278 | 0.126 | 0.1291 | 0.1398 | 0.1407 | 0.1247 | 0.128 | 0.1235 

500 | 0.1308 | 0.131 | 0.1358 | 0.1458 | 0.1458 | 0.1278 | 0.135 | 0.1316 
| | | 


600 1359 | 0.1% | 0.1425 | 0.1530 .1509 | 0.1363 | 0.140 | 0.1396 


700 | 0.1424 | 0.144 


| 0.1492 | 0.1615 | 0.1559 | 0.1473 | 0.150 | 0.1477 
| | | | 
800 1505 | 0.151 | 0.1559 | 0.1701 | 0.1610 | 0.1604 | 0.160 | 0.1558 











II. 
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It is worthy of note that no hysteresis effects of the kind observed by Jaeger 
and Rosenbohm (7) were noticed in this investigation. It is possible that 
these effects occur only when temperatures greater than 800°C. are encountered 
in some part of the temperature cycle. 








Manganese 


The manganese used was in the form of electrolytic nodules, also obtained 
from Johnson, Matthey and Co. The stated amount of impurity was less than 
0.01%, chiefly magnesium. The nodules were slightly discolored by surface 
oxidation. The measured atomic heats are plotted in Fig. 2, and the smoothed 
specific heats are tabulated in Table I, column 5. 
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Fic. 2. Atomic heat of manganese. 





The atomic heat is quite regular from 0° to 700°C., although ithe values are 
considerably greater than 3R. Between 714° and 722°C. a large anomaly in 
thermal capacity was observed. In order to study this anomaly in more detail 
the manganese was heated continuously through this region several times. 
This showed that the anomaly on heating is a sharp phase transition, taking 
place at 717°C. It is attributed to the known transition from the a to the 8 
form. From the continuous heating curves the latent heat of the transition 
is estimated to be 450 cal. per mole. Naylor (11), from measurements of the 
total heat content made by the drop method, estimated the heat of transition 
to be 500 cal. per mole. With continuous cooling, at a rate of about 1° per 
minute, the reverse transformation commenced at about 710°C. and continued 
gradually down to 650°C. It was then completed with an abrupt evolution 
of the remainder of the latent heat, a typical supercooling phenomenon. Al- 
though the a—6 transition in manganese is well known to exist, and has been 
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confirmed by X-ray studies, many different transition temperatures have been 
reported, ranging from 680° to 820°C. The results seem to depend critically 
on the purity and state of the metal. 


The only previous measurements of specific heat of any degree of accuracy 
over the range 0° to 800°C. are those of Wiist ef al. (17). As in the case of 
chromium, no specifications of purity were given. These results are listed in 
Table I, column 6, and the agreement with the present measurements is only 
fair. Wiist et al. did not report any transition between 700° and 800°C. 
Naylor’s measurements of mean heat capacities were not sufficiently accurate 
to take differences over small temperature intervals. 


In measurements taken over a limited temperature range, Ashworth (2) 
observed an anomaly at about 375°C., of the type usually associated with 
Curie point transitions. In order to look for this effect a continuous heating 
curve was taken over the interval 340° to 395°C., but no evidence of any tran- 
sition was found. It is tentatively suggested that Ashworth’s observation 
may have been due to the presence of an impurity which formed a ferromagnetic 
or disordered intermetallic compound with the manganese. 


Cobalt 


The cobalt used was Hilger spectroscopic metal, prepared electrolytically, 
and having less than 0.01% impurity. The atomic heats are plotted in Fig. 3, 
and smoothed specific heats are given in Table I, column 7. The atomic heat 
is quite regular except over the interval 447° to 478°C., where the structure of 
the metal changes from a hexagonal close-packed to a face-centered cubic 
structure. Two continuous heating curves were taken through the region of 
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the specific heat anomaly, in order io determine the shape of the transition 
curve. Atomic heats derived from the continuous heating are shown in the 
inset in Fig. 3. By integrating the peak in this curve the heat of transition 
is estimated to be approximately 60 cal. per mole. 


Cooling curves were taken from 480°C. to considerably lower temperatures, 
but the reverse transition was not detected. This is consistent with other 
studies of the transition. Edwards and Lipson (6) found that cobalt upon 
cooling retains the cubic structure to about 350°C., and then its structure 
gradually changes to the hexagonal form as the temperature is further reduced, 
until about 80°% has been transformed at room temperature. The evolution of 
the heat of transition is therefore spread over a large range of temperature, and 
could not be observed with the limited accuracy obtainable in cooling curves. 


Comparison with other measurements of the specific heats of cobalt is of 
dubious value, since any previous determinations have been made by the drop 
method, with metal of doubtful purity. However, measurements made by 
Wiist e¢ al. (17) and by Umino (16) have been included in Table I, columns 
8 and 9. 


Atomic Heats of the Chromium—Nickel Group of Metals 


The atomic heats of the transition metals, from chromium to nickel, are 
interesting in that they are all abnormally large at high temperatures. In the 
case of the ferromagnetic metals, a portion of the atomic heat can be attri- 
buted to the ferromagnetic energy. This effect disappears at temperatures 
above the Curie point of the metal. 


For theoretical discussion the measured values of the atomic heat at constant 
pressure, C,, have to be corrected to constant volume by the well known 


C,— Cy= &VT/Jx, 


formula 


where a is the coefficient of volume expansion, V the atomic volume, and «x the 
compressibility. For most metals this correction can only be estimated, on 
account of the lack of elastic data at high temperatures. This has been done, 
using figures given by Stoner (14) for nickel, and the best available data 
for a and « of the other metals. In Table II are given what are believed to 
be the most likely values of C, for the five metals concerned. For mangan- 
ese, chromium, and cobalt these have been derived from the measurements 
described in this paper. The values for iron have been derived from meas- 
urements by Wiist e¢ al. (17); values for nickel from those of Sykes and Wil- 
kinson (15). In all cases the corrected atomic heats are still much larger than 
the classical C,= 3R at the higher temperatures. 


According to Debye’s theory of specific heats (12), the portion of the atomic 
heat which is due to vibrations of the crystal lattice should increase from 0 at 
the absolute zero of temperature to 3R at high temperatures. The shape of 
the atomic heat - temperature curve is similar for all solids, the rate of increase 
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depending upon a single parameter, the Debye characteristic temperature 90. 
It is known that the Debye specific heat formula is only an approximation, but 
it is expected to be sufficiently accurate for extrapolation purposes at temper- 
atures JT ~ 9. Values of 6 derived from specific heat measurements at lower 
temperatures are available for manganese (1), iron (9), cobalt (5), and nickel 
(8). Using these and an estimated value for chromium (10) the contribution, 
C,, of the lattice vibrations to the atomic heats has been calculated by the 
Debye approximation. Subtracting the estimated lattice contribution from 
the experimental atomic heats leaves the residues, C, — C,, tabulated in Table 
III. These residues must represent some other form of internal energy, and 
an explanation must now be sought. i 












TABLE II 










ATOMIC HEATS AT CONSTANT VOLUME, Cy (CAL. PER MOLE PER °C.) 











l era 
Temperature, ies, Chromium) Manganese Iron Cobalt | Nickel 


0 | 5.38 | 5.95 | 5.81 5.67 | 6.12 
















100 5.93 6.39 6.47 
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In the case of a metallic crystal there is a contribution to the energy from the 
conduction electrons. It can be shown (10) that at low temperatures this 







contributes an amount 
1 9LTo7Tr/s , , 
to = ns wk?I {/& = al 






per electron to the atomic heat. ) is very closely equal to the highest energy 
level of the conduction electrons which is filled at absolute zero, and the 
approximation is valid when kT << e,. As the temperature approaches the 
value T)= «,/k, known as the degeneracy temperature, the relation is more 
complicated. At high temperatures the electronic specific heat approaches the 








. . . . 3 . . 
classical limiting value = qgR, where gq is the number of conduction electrons 


- 





per atom. The value of g is not known precisely for any of the transition 
metals. However, for approximate calculations it can be taken equal to the 
number of “holes’’ in the d-band of electronic energy levels, because these 
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“holes” are responsible for the greater part of the specific heat, and the con- 
ducting electrons in the s-band of levels contribute only at very high tem- 
peratures. 


The coefficient a in the electronic atomic heat has been determined for nickel 
from low temperature specific heat measurements (8), with the value 17 X 1074 
cal. per degree. From this the degeneracy temperature, T,, can be deduced, 
leading to the unusually low value of 3000°K. Using this value for To, the 
electronic specific heat has been estimated as a function of temperature from 
the curve given by Mott and Jones (10), giving the values in the last column 
of Table III. It appears, therefore, that the high values of the atomic heat 
of nickel above the Curie point can be explained by the effect of the conduction 








electrons. 
TABLE III 
RESIDUES OF THE ATOMIC HEAT, Cute AFTER SUBTRACTING THE 
CONTRIBUTION OF THE LATTICE VIBRATIONS 
Temperature | Chromium] Manganese} Iron Cobalt Nickel | eee ar 
kath 0 eo 4 0.22 * 0.56 ; 0.60 ae 0.43 0.78. ew 
100 0.49 0.75 0.89 0.50 0.86 
200 0.63 | 1.04 1.27 | 0.66 | 1.29 
300 | 0.78 1.13 | L.t¢ 0.94 2.00 
400 0.79 1.31 2.22 | 1.10 1.23 1.09 
500 | 0.91 1.53 | 2.87 1.18 | 1.20 “ 
600 1.12 1.85 3.44 1.55 1.23 1.13 
700 1.44 2.27 4.99 2.15 | 1.24 - 
800 1.82 2.58 2.40 | 2.85 1.26 | 1.16 
Curie cot BOR: + aoc een ie ie ea ea 
temperature | - - 770° —1100° 358° 

















Cobalt and iron both display large residual atomic heats, the greater portion 
of which can be attributed to the internal energy associated with the ferro- 
magnetic spontaneous magnetization. For cobalt the ferromagnetic contri- 
bution.is small below 400°C., and the observed excesses at these temperatures 
are consistent with an electronic term a = 12 X 10, as found by Duyckaerts 
(5). No further conclusions can be drawn for these metals. 


- The electronic term in the specific heat increases linearly with temperature 
at low temperatures, and then at a slower rate as JT becomes comparable with 
Ty. The residual atomic heats for chromium and manganese, as given in Table 
III, exhibit quite a different behavior. The rate of increase at the lower 
temperatures is consistent with an electronic origin, but at the higher temper- 
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atures the values increase more and more rapidly. For manganese the linear 
electron term has been estimated from low temperature measurements to be 
approximately 0.0035 cal. per degree, but the degeneracy temperature cannot 
be satisfactorily estimated because the number of electrons is not known. If 
the degeneracy temperature is of the order of 3000°K., as for nickel, the excess 
specific heat at room temperature should be about twice the observed value. 
If, as seems likely, the degeneracy temperature is very much lower, the high 
excess specific heats at high temperatures are quite unexplained. It appears, 
therefore, that there must be some other source of internal energy besides the 
conduction electrons, and that there are at least two factors contributing to 
the observed results. Since chromium and manganese exhibit similar be- 
havior, it is likely that the same explanation will serve the two cases. 


‘ 


A possible explanation lies in the phenomenon of ‘“‘antiferromagnetism’’ 
which has been observed in certain manganese compounds (4, 13), and which 
has been suggested as a possible explanation of the paramagnetic behavior of 
chromium and manganese metals (3). In these cases the exchange interactions 
cause a co-operative antiparallel alignment of the electron spins in the crystal, 
in contrast with the parallel alignment which gives rise to ferromagnetism. 
The antiparallel orientation of the spins at lower temperatures gives rise to a 
specific heat anomaly of the same type as that shown around the Curie point 
of ferromagnetic materials. Therefore, if chromium and manganese have anti- 
ferromagnetic Curie points at some temperature above 800°C., the excess 
specific heats below 800°C. would be-as they have been observed in this 


investigation. 


Chromium is not known to exhibit any anomalies below its melting point, 
but it may be significant that there is an unusually large decrease in the atomic 
heat at the melting point (16). Naylor (11) has observed a complex series of 
structural changes in manganese between 1000° and 1100°C. However, his 
measurements indicated little change in the atomic heat. In conclusion, it 
should be mentioned that no such anomaly has yet been observed for a pure 
metal, and that this explanation for the abnormal rise of specific heat in 
manganese and chromium is completely unverified. The question cannot be 
settled until magnetic measurements are made in the temperature region 
concerned. 
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PREPARATION AND CHARACTERISTICS 
OF A POLONIUM-BERYLLIUM NEUTRON SOURCE! 






By J. W. T. Spinks? anp G. A. R. GRAHAM? 






Abstract 






The preparation of a polonium—beryllium sandwich containing 2.1 c. of polo- 
nium is described. The source emits 2 X 10° neutrons per second and its Y 
radiation is approximately equivalent to that of 0.2 mgm. radium in equilibrium 
with its decay products. Polar diagrams for the angular distribution of the fast 
neutrons are given. Emission data for two similar neutron sources, containing 
3.2 and 3.5 c. polonium respectively, are given. 












Introduction 


Polonium—beryllium neutron sources are usually made by the process of 
adding a polonium solution to beryllium powder, drying the mixture, and com- 
pressing it into a pellet. This is an unpleasant operation at best and becomes 
somewhat hazardous when large quantities of polonium are handled (owing to 
the ‘“‘volatility’’ of the polonium). The health hazard has had the effect of 
limiting the amount of polonium used in a single source and up until a few 
years ago a 100 mc. polonium—beryllium neutron source was considered ex- 










tremely strong. 





The necessity of treating beryllium with polonium solutions can be avoided 
by making a so-called sandwich in which a platinum foil, on which polonium 
is deposited, is sandwiched between two slices of beryllium metal. This sand- 
wich has the disadvantage that only 50%* of the alpha particles expend their 
energy in the beryllium, but has the advantage that it can be made without 
risk of spreading contamination. This procedure has apparently been used on 
previous occasions (see, for example (3)). However, as experimental details are 
not given in the literature, it was thought worthwhile to record our experience 
in making three polonium—beryllium neutron sources having polonium acti- 
vities of 2.1, 3.2, and 3.5 c., respectively (measured calorimetrically, February 
15, 1945). 


The preparation of the polonium—beryllium source containing 2.1 c. polo- 














nium will first be described. 





Materials 


The polonium was on the lower 20 mm. of a platinum foil 5 by 35 mm. by 
0.12 mm. thick. The foil was hung from the lid of a platinum capsule. The 
lid was friction tight and the capsule was sealed, in helium, in a glass tube. 
The beryllium was in the form of two semicylinders, length 27 mm., diameter 






1 Manuscript received in original form June 22, 1949, and, as revised, October 11, 1949. 
Work done in March 1945 in the Montreal Laboratory of the Canadian Atomic Energy 
Project. Issued as N.R.C. No. 2087. 
2 University of Saskatchewan, Saskatoon, Canada. 
3 Department of Physics, University of Rochester, Rochester, New York. 
* With very thin foil this factor might be increased somewhat. 
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10mm. The flat face of one cylinder had a slot in it, 6 mm. wide and 0.15 mm. 
deep, to take the platinum foil. The beryllium cylinder slid into a brass 
capsule and had a thin beryllium disk at each end. The capsule could be closed 
by a screw and then soldered down (Fig. 1). 


(Cc) 


UZ £L2 





(A) 


Fic. 1. Polonium-beryllium neutron source: (A) assembled, (B) semicylinder beryllium with 
slot for platinum foil, (C) plan. 


Method 


In order to avoid contamination, the operation was carried out in a so-called 
“dry box’’ fitted with rubber gloves and a glass lid, the dry box being placed 
in a hood with a good draft and a good light. Before the operation was 
started all necessary tools were placed in the box. The lid of the box was 
then sealed down and nitrogen streamed through for some minutes. The glass 
tube containing the platinum capsule was then broken (hot wire) and the 
capsule removed. Next the slices of the sandwich were slid out of their con- 
taining case, the foil laid between them, and the parts of the cylinder slid into 
the brass capsule. The polonium could be seen as a brownish black deposit 
on the lower part of the foil. The top part of the platinum and lid were snipped 
off and the small overlapping piece of platinum tucked in. The top beryllium 
disk was put in place and the lid screwed in. The capsule and pieces of glass 
tube were placed in a screw top bottle for later disposal. At this stage a test 
with a portable a-counter showed no contamination outside the box. The 
nitrogen was then turned off and the box opened. The capsule was removed 


and sealed down with soft solder. 


The platinum capsule from the 2.1 c. of polonium was shown later to contain 
approximately 74 mc. polonium. This means that about 4% of the polonium had 
‘wandered’ from the foil to the surrounding capsule in the course of a few weeks. 
The gloves and tools which had touched the platinum capsule showed some 
a-activity but the fingers of the operator were not contaminated; in fact, the 
contamination was limited to the inside of the dry box. The whole operation 
took about 15 min. (naturally it had been practiced several times with a 


dummy). 
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Radiations from Polonium-—Beryllium . 
-Radiation.—The y-rays, measured through 1 cm. lead, were equivalent to 
those from 0.2 mgm. radium in equilibrium with its decay products. 





Angular Distribution of Fast Neutrons 


In order to measure the angular distribution of fast neutrons, the source was 
mounted on a turntable, 18.5 cm. from the face of a truncated paraffin cone 
in which was enclosed a boron trifluoride chamber (Type E (5)) (see Fig. 2). 


NEUTRON 
1ccm. SOURCE 


BF CHAMBER 


PARAFFIN BLOCK 








Fic. 2. 
source. 


Arrangement for measuring neutron distribution about polonium—beryllium neutron 





The boron trifluoride chamber was connected to the usual discriminator and 
counting circuit and readings were taken every 30 degrees, first with the source 
mounted with its axis at right angles to the surface of the turntable and next 
with its axis parallel to the surface of the turntable. Approximately 750 counts 
were recorded at each position (probable error 2.7%). 


The results are recorded in Tables I and II and in polar diagrams in Fig. 3. 
It appears that the angular distribution about the axis of the cylinder is 
symmetric but the distribution in a plane containing the axis of the cylinder 
TABLE I 


NEUTRON DISTRIBUTION ABOUT THE CYLINDER AXIS 
(See Fig. 2) 











Angle, degrees* | Counts per 10 minutes 
0 725 
30 718 | 
60 744 
90) | 705 
120 | 796 
150 | 797 
180 725 
270 703 
Probable error for 750 counts is 2.7% 
* Cylinder rotated about the cylinder axis and 
angle measured with reference to a line at right 


angles to the cylinder axis. 
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is slightly asymmetric. The neutron intensity at 20 cm. from the source is 
about 20% greater in the broadside-on position than in the end-on position. 
Such a result is not unexpected for a source of this type. 


TABLE II 


NEUTRON DISTRIBUTION IN A PLANE CONTAINING 
THE AXIS OF THE CYLINDER 








Angle, degrees* | Counts per 10 minutes 
0 (broadside-on)| 816 
| 774 
60 686 
90 (end-on) | 600 
120 | 694 
150 774 
180 (broadside-on ) 77 
270 (end-on) 704 





* Cylinder rotated about the normal through 
the mid-point of the cylinder axis and angle mea- 
sured with reference to a line at right angles to 
both the cylinder axis and the axis of rotation. 


(A) (8) 


Fic. 3. Angular distribution of fast neutrons about an axis: (A) through the cylinder axis, (B) 
at right angles to the cylinder axis. 


Neutron Emission (Neutrons per Second) 


The total neutron emission per second was measured by comparison with 
a 110 mgm. Ra—a-Be neutron source, stated to be equivalent to 64 mc. radon- 
beryllium. The method involved suspending the source and a boron tri- 
fluoride chamber in a tank of distilled water and taking counts for various 
distances between source and chamber (see, for example (1, 2)).. About 2500 
counts were recorded at each position (probable error, 1.3%) and the results 
plotted as density versus radial distance squared, in arbitrary units (Tables 

III and IV and Fig. 4). 


Points on the curve for values of R less than 4 cm. were calculated assuming 
a parabolic relation between JR? and R; points on the curve for large values of 
R were calculated assuming an exponential relation between J and R (e.g., 
relaxation length for the polonium—beryllium source = 6.1 cm.). 
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TABLE III 


VARIATION OF NEUTRON INTENSITY WITH DISTANCE FROM 
A POLONIUM-BERYLLIUM NEUTRON SOURCE IN WATER 













| 


Besa. i Uinta 
Distance R, Intensity J, Intensity J’, IR? 
















cm. counts/min. counts/min. 
4.6 1694 912 | 35800 
5.0 1618 874 40450 
6.0 1405 756 50500 
7.0 1240 670 | 60800 
8.0 1002 540 64100 
9.0 850 459 68800 
1.6 1706 check 
10.0 695 375 69500 
11.0 581 70300 
12.0 482 | 69500 
14.0 312 61100 
17.0 177 51100 
20.0 100 | 40000 
23.0 54 | 28600 
27.0 26 18950 
1.6 1700 check 
36.3 4.5 5920 
4.6 1682 check 







TABLE IV 


VARIATION OF NEUTRON INTENSITY WITH DISTANCE FROM 
RADIUM-BERYLLIUM NEUTRON SOURCE (110 MGM.) IN WATER 





















Distance R, cm. Intensity J, counts/min. IR? 


18430 















4.2 1044 | 

5.0 930 23260 
6.0 834 | 30000 
7.0 670 | 32800 
8.0 540 34600 
9.0 438 | 35500 
10.0 345 | 34500 
11.0 280 | 33900 
12.0 224 32300 
4.2 1058 check 
4.2 1060 check 
14.0 144 28200 
17.0 78 | 23100 
20.0 40 | 16000 
23.0 22.8 12090 
27.0 | 11.2 8150 
4.2 1070 | check 
4.6 cm. 1682 check 

(Po-Be) 
Area under curve for Po—Be = 129.8 * 104 
Area under curve for Ra—a—Be = 62.2 X 10! 
ae : : 129.8 F 
Thus the Po—Be source is equivalent to me xX 110 = 230 mgm. Ra-a—Be 
Va. 


9. 
or to < 8 xX 64 = 134 mc. Rn-Be. 
12.2 
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UPPER CURVE Po-Be 
LOWER CURVE Ra-Be 
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Fic. 4. Distribution of neutron density I multiplied by R?® for polonium-beryllium and 
Ra-a-Be neutron sources in water. 


Assuming that 1 mc. Rn—Be emits 15,000 neutrons/sec., the 


Po-Be source emits 2.0 X 10° neutrons/sec. 


(March Ist, 1945.) 

Two further polonium—beryllium neutron sources were prepared by the same 
method. Their characteristics, together with those of the No. 1 source, are 
summarized in Table V. 

TABLE V 


CHARACTERISTICS OF THREE POLONIUM—BERYLLIUM NEUTRON SOURCES 
(as oF Marcia 1, 1945) 


| Ra-y-ray equivalent 








Number Curies Po Neutrons/sec. |(through 1 cm. lead), mgm. 
ry 1 . 2.1 2.0 X 10° 0.2 
2 3.2 3.0 X 10° 0.3 
3 3.5 3.2 X 108 0.35 


In Fig. 5, density is plotted versus distance for points close to the source. 
The actual values for the radium—beryllium neutron source have been plotted, 
while the values for the polonium—beryllium neutron source have been nor- 
malized to the same intensity (540 counts per min.) at R = 8.0 cm. (J’ in 
Table III). The greater slope for the radium—beryllium source indicates a 
greater abundance of low energy neutrons for it than for the polonium-bery]- 


lium source. 
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Fic. 5. Plot of neutron density versus distance, close to source (immersed in water). 


From a plot of /R* versus R, the area [ear may be obtained. Dividing 


this by the area under the curve for 7R? versus R allows the calculation of 
‘ry? average ’; 


| 
sr| IR‘dR 
yi = = ; 
ir TR°dR j 
0 

The migration length is defined by Ly = / r . For the polonium—beryllium 
6 ‘ 

source, 7? = 351 cm.? and Lo= 7.6 cm. Seligman (4) gives 7? = 330 cm.? and 


Lo= 7.4 cm. for radium—beryllium and a boron trifluoride detector. 
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A MAGNETIC FIELD STABILIZATION CIRCUIT! 


By L. Katz, P. A. Forsytu, L. F. CupNEy, G. W. WILLIAMS, 
H. E. JoHNs, AND R. N. H. HASLAM 


Abstract 


The paper describes the apparatus used to stabilize the magnetic field of a large 
electromagnet. A rotating coil placed in the magnetic field generates an alter- 
nating voltage which is proportional to the strength of the field. Variations in 
this rotor voltage are used to control the current supplied to the electromagnet 
in such a manner as to keep the field constant. By means of a calibrated field 
control the magnetic field strength may be adjusted to any value between 500 
and 6000 gauss. The current for the electromagnet is obtained by rectifying 
60 cycle alternating current in a full wave, thyratron, rectifier circuit. This 
current is controlled by a pulse circuit which governs the fraction of each cycle 
during which the thyratrons are conducting. The stabilization at any field 
setting is better than 1 part in 1000. 


The auxiliary apparatus under development for use with the University 
of Saskatchewan’s new 22 Mev. betatron includes a beta-ray spectrometer 
which will be used to measure the energy distribution of high energy beta 
particles. The spectrometer is of the plane field, 7 radian type in which the 
particles are bent through 180° by the action of a strong magnetic field. In 
order that precise measurements of electron energies may be made, it is neces- 
sary that the magnetic field of the spectrometer be maintained at a constant 
value over relatively long periods of time. Accordingly the construction of a 
suitable field stabilizing circuit was undertaken early in 1948. 


The problem which presented itself may be summarized as follows: 


(i) A power supply was required which could provide a maximum of 1.5 
amp. at 600 v., to the electromagnet. 


(ii) The field of the electromagnet was to be maintained at a constant value 
within approximately + 0.1% over a period of several hours. 

(iii) The field was to be continuously variable over a range of at least 10 
to 1, by means of a simple control. 


(iv) The control of (iii) was to be calibrated in terms of field strength. 


It was decided that (i) could be satisfied by the use of a full wave rectifier 
circuit, using the large inductance of the electromagnet as part of the filtering 
circuit. In the final instrument, thyratrons are used as the rectifying tubes. 
Control of the output current is obtained by changing the a-c. voltage, and 
regulation is obtained by varying the portion of the cycle over which the 
thyratrons are conducting. This method of control provides a distinct saving 
in power over the more conventional series tube regulator. 

1 Manuscript received August 19, 1949. 


Contribution from the Department of Physics, University of Saskatchewan, Saskatoon, 
Saskatchewan. 
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In order to fulfill the requirements of (ii) and (iii) and particularly that of 
(iv), it was decided to build a stabilizing circuit which was sensitive to the 
instantaneous field strength rather than one which provided current stabili- 
zation. This method of stabilization avoids the necessity of compensating for 
the rather large hysteresis of the iron core. 


In the apparatus which has been constructed a rotating coil is placed in the 
air gap of the electromagnet and driven by a synchronous motor. The ampli- 
tude of the alternating voltage produced by this rotating coil is taken as a 
measure of the field strength. A fraction of this voltage is balanced against 
the voltage produced by a similar coil rotating between the poles of a strong 
permanent magnet. Any variation in the field strength of the electromagnet 






REFERENCE VOLTAGE 






PERMANENT 
MAGNET 







THYRITE PROTECTIVE 
RESISTORS 


SPECTROMETER 
ELECTROMAGNET 






Fic. 1. Schematic diagram. 


destroys the balance and results in an error signal which is transmitted to the 
stabilizing circuit. The stabilizing circuit is designed to vary the current sup- 
plied to the electromagnet in such a way as to reduce this error signal to zero. 
Control of the magnitude of the field is then provided by varying the fraction 
of the electromagnet rotor voltage which is balanced against that of the perma- 
nent magnet rotor. 


A schematic representation of the apparatus is given in Fig. 1, from which 
the stabilization loop can readily be traced. A simplified circuit diagram is 
given in Fig. 2. The 30 cycle alternating voltages produced by the two rotors 
are connected out of phase to the two ends of a potentiometer. The error 
signal is produced between the slider of the potentiometer and ground, and is 
applied to the input of the 30 cycle tuned amplifier. This amplifier consists 
of three stages and has an over-all gain of 300. The tuning is accomplished 
by means of a twin-T (4, 5) connected between anode and grid of the first 
stage, V;. The last stage is push—pull in order to provide paraphase voltages 
for the next unit of the circuit. 
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The amplified error signal is detected by means of the phase sensitive de- 
tector. This circuit produces a d-c. voltage, the amplitude being determined 
by the degree of unbalance between the two rotor voltages, and the polarity 
indicating which is the greater. The fact that this circuit is balanced with 
respect to ground ensures a very stable rest point which is independent of the 
amplitude of the reference signal. 


The integrator is similar in design to the conventional ‘‘memory”’ circuits 
which have been treated in the Radiation Laboratory Series (3). Its purpose 
is to integrate the output of the phase sensitive detector. In the presence of 
an error signal the output potential of the integrator drifts in either the posi- 
tive or negative direction, depending on the phase of the error signal. As the 
error signal is reduced to zero it approaches a new stationary value. The inte- 
grating condenser is connected between the cathode of V3 and the control grid 
of V;, while the integrating resistance is provided by the internal resistance of 
the phase sensitive detector. This arrangement gives an effective time con- 
stant of five seconds. 


The remainder of the stabilization loop is designed to vary the point of firing 
of the thyratron rectifier tubes in accordance with the potential which has 
been established at the output of the integrator circuit. When this potential 
is increased the thyratrons are caused to fire later in the cycle and hence the 
current through the electromagnet is reduced. 


The trigger pulse circuit which consists of the overdriven amplifiers, 12 and 
V13, and the associated diodes, produces a sharp pulse each time the 60 cycle 
supply voltage passes through zero, 1.e., 120 pulses per second. These pulses 
are used to trigger the phantastron delay circuit which produces another series 
of pulses delayed with respect to the trigger pulses by a time dependent on 
the potential existing at the cathode of Vs. This type of circuit has been 
treated extensively in the literature (see for example Close and Lebenbaum (2) 
or B. Chance ef a/. (1) ). In addition to the potential at the cathode of Vx, the 
phantastron delay time is sensitive to the potential to which the control grid 
resistor of Vi, is returned. This method of control is used in conjunction with 
the voltage feedback loop. 


The delayed pulses from the phantastron are amplified in the pulse amplifier 
and applied to the grids of the thyratron rectifier tubes. A self-biasing ar- 
rangement is provided in order to prevent the thyratrons from firing in the 
absence of the triggering pulses. The gas regulator tube in this bias circuit 
ensures that the bias never becomes greater than the amplitude of the trig- 
gering pulses. In adjusting the apparatus to provide a certain field strength 
the a-c. voltage applied to the thyratrons is set by means of the variac control 
to a value which supplies approximately the correct field without the stabili- 
zation action. This procedure ensures that the corrections required of the 
stabilizing circuit will be small and the thyratrons will always be firing in 
approximately the same portion of the cycle. 
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The condenser connected across the terminals of the rectifier operates in 
conjunction with the large inductance of the electromagnet (approximately 
250 h.) to smooth the current passing through the magnet coils. A thyrite 
resistor is connected across each of the magnet coils in order to provide pro- 
tection against overvoltage in event of a power failure. 


The rate at which the field is returned to its correct value by the stabilizing 
circuit, as described above, is proportional to the amplitude of the error signal. 
In order to prevent oscillations it is necessary to introduce a second stabili- 
zation loop which is sensitive to the rate of change of the field. Fortunately, 
owing to the inductance of the electromagnet, the a-c. component of the poten- 
tial drop existing across its terminals is proportional to the rate of change of 
the field. Half this voltage is passed through the voltage feedback amplifier, 
Vis, to the phantastron. In this way, long time corrections are applied to the 
phantastron by the integrator circuit while rapid fluctuations in line voltage 
are compensated for by the voltage feedback circuit.. The time constant of the 
Vis control grid circuit is made sufficiently long to ensure a suitable overlap 
between the two correcting actions. This grid is also shunted by a condenser 
in order to reduce the gain of the circuit at a frequency of 120 cycles since a 
fairly large ripple of this frequency is present in the output of the rectifier 
circuit. 

The switching arrangement shown between the cathodes of Vs and V4 pro- 
vides ‘‘start’’ and ‘‘run” positions. In the ‘‘start’’ position the cathode of V» 
is disconnected from the integrator circuit and connected to a fixed potential 
chosen so that the thyratrons will be caused to fire near the peak of the supply 
voltage. The zero-center meter in the circuit indicates the difference between 
the integrator output potential and the fixed voltage on the cathode of V5. 
This arrangement allows the unstabilized value of the field to be adjusted 
initially, by means of the variac control, so that only a small transient will 
result when the stabilization loop is completed by switching to ‘‘run.”’ In the 
latter position a large positive or negative meter reading indicates that the 
circuit is operating near the limit of stabilization. Another section of the same 
switch disconnects the voltage feedback amplifier in the ‘‘start’’ position in 
order that it will not oppose adjustments to the field strength made by means 
of the variac control. 


In addition to the thyrite protection resistors a number of interlocking 
switches and relays, not shown in Fig. 2, have been incorporated into the 
apparatus. These serve to prevent damage to the magnet or associated cir- 
cuits as a result of power failure or inadvertent operation of the controls. 
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THE ELECTRON CYCLOTRON! 
By P. A. REDHEAD, H. LECAINE, AND W. J. HENDERSON 


Abstract 


A magnetic resonance accelerator for electrons is described which was con- 
structed on the principles suggested by Veksler. A constant magnetic field and 
applied frequency are employed, the frequency of the accelerating radio-frequency 
field being 2800 megacycles per second. Final energies of five million electron 
volts are obtained. 


1. Introduction 


In a resonance accelerator such as the cyclotron, the rotational frequency 
of the accelerated particles is matched to the frequency of the applied radio- 
frequency field. The relativistic increase in mass of the accelerated particle 
as it gains in energy causes a change in its rotational frequency, and thus the 
particle falls out of synchronism with the radio-frequency field. This process 
imposes an upper limit on the particle energy obtainable from this type of 
accelerator. This effect occurs at thirty to forty million electron volts for 
heavy particles such as protons, but at quite low energies for electrons. For 
this reason it has never been possible to apply the conventional cyclotron 
technique to the acceleration of electrons. 


To overcome this relativistic limitation, Veksler (7) and McMillan (5) inde- 
pendently proposed the principle of phase stable orbits in magnetic resonance 
accelerators. Veksler proposed two types of accelerators. One of these pro- 
posals involves the concept of phase stable orbits in which the particles gain 
energy as a result of a slow variation of magnetic field, or the applied radio- 
frequency, or of both. Accelerators of this type are the “Synchrotron” arid 
the ‘‘Synchrocyclotron’’. The second of Veksler’s proposals utilizes a constant 
magnetic field and applied frequency, and is the subject of this paper. Since 
this accelerator resembles the cyclotron in many respects we have called it 
the “Electron Cyclotron”’. 


The period of rotation of a particle in a magnetic field is a linear function 
of the total energy of the particle. Thus, when the particle is periodically 
accelerated, the difference in times taken to perform successive orbits is a 
constant, independent of the total energy, and is a function only of the kinetic 
energy gained by the particle in crossing the accelerating gap. Resonance 
acceleration can therefore be obtained if this constant time difference is made 
equal to an integral number of periods of the accelerating field. In this case 
the particle will always arrive at the accelerating gap at the same phase of the 
radio-frequency field. An initial condition is required to ensure that the time 
of revolution of the particle in the first orbit is an integral number of periods 

1Manuscript received August 8, 1949. 


Contribution from the Radio and Electrical Engineering Division and the Physics Division, 
National Research Laboratories, Ottawa, Canada. Issued as N.R.C. No. 2058. 
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of the accelerating field. The resonance conditions are derived in the next 
section. 

An accelerator has been constructed on these principles and acceleration of 
electrons to 4.6 million volts has been achieved (2). Fig. 1 is a simplified plan 
view of the accelerator through the median plane. Electrons produced by field 
emission from the inside surface of the resonant cavity are accelerated across 
the gap in the cavity to a maximum of about 600 kv. These electrons emerge 
from the hole in the pole of the cavity and rotate in the magnetic field of 
approximately 1000 oersteds. The electrons returning to the accelerating gap 
in the correct phase are once more accelerated, gaining about 600 kv. at each 
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Fic. 1. Simplified plan view of accelerator through the median plane. ‘‘A” is an alternative 
target. 


acceleration. Since the energy gained at each crossing of the gap is so large, 
the orbits, which are cotangential at the accelerating gap, are widely spaced 
from one another at the outer portion of the orbits. Eight orbits are possible 
in the present 14-in. diameter vacuum chamber. 

The cavity, which is ellipsoid-hyperboloid in shape, is excited from a 
tunable magnetron oscillating at 2800 Mc. per sec. The magnetron is operated 
at 2u sec. pulse-width, and a repetition frequency of 435 c.p.s. The radio- 
frequency system is described in more detail in Section 8. 


2. Relativistic Equations of Motion 


For circular motion in an axially symmetric magnetic field, with an external 
electric torque applied about the axis of symmetry, we find that cylindrical 
co-ordinates give the simplest form of the equations of motion. Neglecting 
second order terms and using rationalized m.k.s. units throughout, these equa- 
tions may be written as: 


oe .. mr& — er6B, (la) 
dt 

d(ms) _ «iB, (18) 
dt 

oe = —erB + Ty, (1c) 
at 


where m is the relativistic mass of the electron, 


“ 
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B, and B, are the axial and radial components of the magnetic field at 
the orbit, and 


T, is the applied external torque about the axis of symmetry of the orbit. 
It is given by: 
qT, = dW dL (2) 


=> —— 


dé dé 


dw . ; ins . 
where —— is the applied electrical torque, dW being the work done by this 
torque in an angular displacement dé, and 


dL. 4 aie 
= is the decelerating torque due to radiation from the electron. 
t 


The radiation loss (3) depends on the radial acceleration of the electron and 
is given by the following relation: 


2(77)4 
at = Se {U\", (3) 
3R (Uo 
where AL is the energy loss per revolution, 
U is the total energy of the electron, and 
Us is the rest-mass energy of the electron (511.24 X 10° electron volts). 


This term is significant only for electrons with energies greater than several 
hundred Mev. 


The conditions for a stable orbit are: 
r=R;2=0;F =0; B,=0; 6=0; B,= B. (4) 


Then, from (1a) 
ewRB = mRo 


QO = — = - (5 


To find the orbit radius for a given electron energy and magnetic field, we make 
use of the well known relativistic relation: 


pm (me) wo SO Se) (6) 


Substituting in Equation (5) we obtain 


v _ (U?— Us)? 


Ro ese 


w eBc 


(7) 





This equation may be rearranged to give a convenient expression for the 
electron energy (V) in electron volts, in terms of the orbit radius and magnetic 


field. 


V = 300 4/ (BR) 4. (moc) _ 5.1124 x 105 (8) 
€ 
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9.107 XK 10~% gm. 
2.9978 X 10!° cm. per sec. 
e = 4.802 X 107 e.s.u. 


ll 


B in oersteds Mo 
R in cm. 


o 
| 


3. Conditions for Resonance 


If the time of revolution of an electron in the uth orbit is T,,, then 


ry = 28 2 2Us 


w eBc? 





from Equation (5). (9) 


If the electron is accelerated through a potential drop of V volts at each 
crossing of the accelerating gap, then 
U, = Uot neV. 

_ 24V 
Be? 
Thus it is seen that the result of the relativistic increase of mass is to intro- 

duce a constant difference in the times taken to perform successive orbits. 

This time difference is independent of the total energy of the particle and in 

a constant magnetic field it depends only on the voltage on the accelerating gap. 

This important concept was first presented by Veksler. 





Thus, eee «&f (10) 


For resonance acceleration of the electrons, this difference in the time taken 
to perform successive orbits must be an integral number of periods of the 
applied radio-frequency field. 


Thus, AT = wcll = br, (11) 


2 


CB 
where 0 is any integer except zero, 
7 is the period of the radio-frequency field, and 
V’ and B’ refer to the resonance values of voltage and magnetic field. 


Equation (11) is the first resonance condition relating the time difference 
between successive orbits and the radio-frequency period. A second condition 
is necessary to ensure that the initial entrance energy of the electron is such 
that after acceleration at the gap it will perform the first orbit in an integral 
number of periods. This condition is fulfilled when 

Tin 2n(Uo + eVi + eV") _ ‘ae, (12) 
€B’c? 
where a is any integer except zero or unity, and 
V; is the injection voltage of the electron. 

When I; is zero, and e€V’= U, it will be seen that a = 2b. Thus resonance 
acceleration can be achieved with electrons starting from rest if the accelerating 
potential is made equal to the equivalent rest-mass potential of the electron 
(511.2 kv.), and if the frequency and magnetic field are chosen such that 


ch (13) 
eB’ c? 


| 
j 
| 
| 
| 
| 











REDHEAD ET AL.: THE ELECTRON CYCLOTRON 77 


To avoid the problem of high energy injection this mode of operation was 
decided on, and thus for the accelerator described in this paper Equation (12) 


reduces to on (Us + <V’) 


= ar. 14 
eB’ c? (14) 





Equations (11) and (14) define the resonance conditions, and from these 
equations we find that 
ae _ 20m (15) 

e(a—b)r 


and ve = moc* aes 
é (¢@—5) 


From consideration of available radio-frequency generators, it was decided 
to operate in the a = 2, b = 1 mode where the resonance parameters have the 
following convenient values: 


(16) 


B’= 995 oersteds 
V’ =. SEL ky. 
f = 2800 Mc. per sec. 


4. Motional Equations for an Electron in an Accelerating Gap 


The choice of the length of the gap across which the electrons are accelerated 

is affected by the following considerations— 

1. The peak voltage developed across the gap for a given power input is de- 
termined by the shunt resistance of the cavity. The shunt resistance in 
general decreases with gap length. Thisisdiscussed more fully in Section 9. 

2. The transit time of the electron across the gap will be an appreciable frac- 

tion of the radio-frequency cycle. The transit time, which is proportional 

to the gap length, must be short enough to allow the electron to traverse 
the gap during the positive part of the radio-frequency cycle. 

3. The foregoing analysis of the electron motion has assumed that the elec- 
trons are accelerated across an infinitely short gap. A gap length that is 
too long will cause displacement of the orbits. 


As a compromise amongst the above factors a gap length of 1 cm. was chosen. 
Having chosen the length of the accelerating gap a calculation of the energy 
gained by an electron in crossing this gap can be made. 


If an electron crosses a gap of width d centimeters, across which exists a 


voltage : ‘ 
8 V = Vo cos wol, 


then the momentum of the electron at time ¢, neglecting the effect of the 
magnetic field, is given by 





t 17 
i - | "cos wotdt + pi, (17) 
ty 











78 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


where /; is the momentum of the electron on entering the gap, and ¢; is the time 
at which the electron enters the gap. Then, 


p — — ld eVo {sin C71) pom sin ¢:} + Di, (18) 


yi v2 wod 
A aa 
c 





where ¢ = wot 
di= Wot. 


Thus an approximation for the momentum of the electron on leaving the 
accelerating gap is, 
eVo : 
po = —cos a2 sin x/2 + p; (19) 


0 


where x is the transit angle, 
a is the phase angle of the radio-frequency field when the electron 
crosses the centre of the gap. 
For a constant transit angle, Equation (19) may be rewritten: 


pbo= A cosa + pi, 
where A= me sin x/2. (20) 


Equation (18) may be solved for v and graphically integrated for given 
entrance angles and entrance velocities to obtain the time taken by electrons 
to crossa lcm. gap. Once the transit angle is known, the final exit velocity 
may be found from Equation (18). 


The results of these calculations for a 1 cm. gap and an injection potential 
of 10,000 v. are shown in Figs. 2a and 2b. These curves correspond to different 
values of the peak accelerating voltage. It will be observed from these curves 
that electrons starting with low velocities and entering the gap at the correct 
phase will obtain almost the peak gap voltage. An experimental check of the 
voltage acquired by electrons starting from rest was obtained by measuring 
the energy spectrum of the first orbit (see Fig. 3). This curve was found by 
plotting the current to a movable collector with a small aperture (1 mm. wide) 
in a known magnetic field, against radius. Conversion to an energy scale was 
achieved by use of Equation (8). The maximum voltage of the electrons in 
this spectrum corresponds closely to the estimated peak gap voltage deter- 
mined from measurements of power input and Q of the cavity. 


The fact that a low-velocity electron could obtain such a high percentage of 
the peak gap voltage was unexpected, and further encouraged the construction 
of an accelerator with zero injection energy. 


On later crossings, where the velocity of the electron approaches the velocity 
of light and may be considered constant, the energy gained on crossing the gap 
is the average of the cosine curve between the entrance and exit angles. The 


| 
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transit angle for a hyper-relatavistic electron crossing a 1 cm. gap is 34°, and 
it can be shown that within the phase stable region the energy gained varies 
between 94% and 99.5% of the peak voltage on the gap. It may be seen from 
Fig. 2a that the maximum energy gained by an electron on the first crossing 
for a peak voltage of 500 kv. is approximately 97% of the peak voltage. Thus 
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Fic. 2(aandb). Variation of voltage efficiency and transit angle with entrance angle for differ- 
ent values of the peak voltage on the gap. Curves are plotted for an entrance energy of 10\kv. ¥ Effi- 
ciency is expressed as a percentage of the peak voltage. 
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an electron starting from rest will gain almost the same energy as an electron 
in the outer orbits traveling with the velocity of light. 


12 


CURRENT (ARBITRARY UNITS) 


0 “| 2 3 “4 5 6 7 8 ‘9 rr) 
ENERGY (Mev.) 


Fic. 3. Current versus energy in the first orbit. 


5. The Principle of Phase Stability 


The inherent phase stability of the electron cyclotron can best be understood 
by a qualitative description of the electron motion. 


Consider first the electron which crosses the gap at a phase ¢p slightly after 
the peak of the radio-frequency cycle (see Fig. 4). This electron will gain the 
resonance voltage V’ (neglecting transit-time effects), and if the magnetic field 
is adjusted to the correct resonance value this electron will arrive back at the 
gap at the same phase @» and receive the same voltage V’. This electron is 
called the resonant electron. Next consider an electron crossing the gap ata 





Seite ae! 
_ & me 





' Oo 


Fic. 4. Illustration of phase stability. 


phase ¢; a little early with respect to the resonant electron. This electron will 
gain more energy than the resonant electron and thus it will take longer to 
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perform the next orbit (see Equation (9)). Thus at the next gap crossing this 
electron will cross at a phase ¢2 slightly closer to the resonant phase. After a 
further acceleration the electron will cross the gap at the resonant phase, but 
with an energy excess, so that it will continue its phase shift to the other side 
of the resonant phase. Now the process is reversed and the electron will be 
returned to the resonant phase position. This process represents a phase 
oscillation about the resonant phase ¢o. Thus there is a phase stable region 
about the resonant phase: within this region electrons will remain in a stable 
orbit; outside this region electrons will be lost from synchronism and spiral 
outwards or inwards and finally strike the walls of the vacuum chamber or the 
cavity. 

Electrons crossing the gap at the correct phase but with an incorrect energy 
will also oscillate about the resonant phase. The oscillations will be very 
similar to the phase oscillations but will differ slightly in phase. 


6. Analysis of Phase Stability 


Proceeding from the gap-crossing equations of Section (4) we may make 
an approximate determination of the phase stable region. From Equation (19) 
it can be shown that the momentum of an electron after the mth crossing of the 
accelerating gap is given by 

Pn = A (COS an + COS an—-1 +... . + COS a1), (21) 
where a, and a,_; are the angles at which the electron crosses the center of 
the gap on the mth and (m — 1)th crossing, and it is assumed that the gap 
transit angle is constant for all crossings. This assumption is reasonable since 
the electron reaches a velocity of 0.9 c. at the first crossing. 

From Equation (6) it can be shown that 


U=c(pe?t+ myc2)'/” 


since p = mv 


fs dU cp 4 “ 
Phen — = p =v U = me 
dp l 
a du dp 
Phus Cm pee 22) 
da da 
Combining Equations (22) and (21) 
dU dp» ‘ j dan— . da 
tes Pe es Ady tein a, +sin an-1— — +... + sin —{ , 
dan day ( day day 


where v, and U, are the velocity and total energy of the electron at the mth 


orbit. 


dan—1 dan—2 


(sin GQn—1 + SiN an—e ee 


dan dan—1 


= — Ap, sin Qn + 


+...+ sina, des )} 


dan—1 











82 CANADIAN JOURNAL OF RESEARCH. VOL. 28, SEC. A. 


* da orl dU ew 
= — Ap, sin a, + —— —— 
dain dage D5 1 


a ee nn (23) 


dan Vn-1 





From Equation (9) it can be shown that 





dT, _ 20 
dU, ¢Bc? 
or dant: _ dan 2nUo (24) 


dU, ei: «Bc? 


Since an41= an + wol'n, where the variation of a, with U, is neglected. 
Combining Equations (23) and (24) we find that 
dU, eBc? Uy 


— = — Av, sin an + — 
day TW0 Un—-1 


(25) 





Considering the extremes of the phase stable region, we see that the condition 


for stability is: E . 
: Al C= = AL n—1 OF Aan = Aika 35 


Thus from Equations (24) and (25) the conditions for stability become 














. «Bc? Un 
Av, SIN an = Be J, beac ‘ 
TO Un-1 
. Bed ; Ur 
sis op (26) 
4rvnVo sin x/2 \ ha 
, * . és mS ‘ Be? V'c 
Using the condition for resonance given in Equation (11), : = 7 where 
Tv 


\ is the wave length of the radio-frequency power and assuming v, = c Equation 
(26) reduces to 
; V’ d/x v oe 
sin a, = — penet ee Sy + — ‘ (27) 
Vo 2 sin x/2 l Un—1) 





Assuming that the resonant phase a’ lies in the center of the phase stable region 
and that one extreme of this region is the peak of the radio-frequency wave, 
7 
then a,= 24a and — = cos Aa’ = cos Aa. Thus Equation (27) becomes 
0 





sin Aa = ca S, a nm \ . (28) 
4 sin x/2 \ e-5) 
ee d ‘> : ‘ 
oe d _ 10.7 = 0.112, x = 40° and for the first two orbits the 
we V1 


following approximation for the phase stable region is found. Aa = 10°. 
Thus the phase stable region extends from the crest of the wave for 20°. 
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Resonance acceleration has been observed over a greater range of magnetic 
fields than the above figure for phase stability would indicate. However, a 
change in the magnetic field causes a shift in that part of the first orbit energy 
spectrum that is utilized for further acceleration and also changes V’, the 
increment of energy gained by the resonant electron. The permissible varia- 
tion in magnetic field for which an output beam is obtained is therefore not a 
function only of the phase stability. 


7. Limitations on the Orbit Current 


It may readily be shown that the maximum attainable current in the accel- 
erated electron beam is determined by the radio-frequency power supplied to 
the accelerating cavity. The electronic conductance of the electron stream 
is effectively in parallel with the shunt conductance of the resonator, thus any 
increase in electron current, and therefore an increase in beam conductance, 
will cause an increase in the effective shunt conductance of the cavity. Fora 
fixed power input to the cavity any increase in shunt conductance will cause 
a decrease in the peak voltage across the accelerating gap. If this process is 
continued far enough the gap voltage will drop below the value necessary for 
resonance acceleration. Thus the electron beam current which lowers the peak 
gap voltage to the resonance value is the maximum orbit current which may 
be employed. 

An analysis of the loading of the cavity by the electron stream allows an 
estimation of the maximum current obtainable with the available radio-fre- 
quency power (250 kw. peak power). 


The power in the cavity when it is not loaded by an electron beam is given by 


W, = V0Gsn, (30) 
2 


where Vo is the peak voltage at the accelerating gap and G,, is the shunt con- 
ductance of the cavity defined as 

2 X energy lost/sec. 

Se ee 

voltage? 
where the voltage is defined as the line integral of electric field along the path 
followed by the electrons. We shall assume that the shunt conductance calcu- 
lated for the cavity is located at the accelerating gap. 
When the cavity is loaded by an electron beam the power in the cavity is 


given by r2 
r Z yy ‘ © 
Wr = — (Gen + Gp), (31) 
2 
where Gz is the electron beam conductance given by 


I (32) 


’ 





Gz = - 


L 
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where V, is the peak voltage on the gap when the cavity is loaded and J is an 
equivalent loading current at the voltage V;. The current in the first orbit, 
which orbit is the major cause of the loading effect, is widely distributed in 
energy. Thus to find the equivalent loading current at the voltage V it is 
necessary to integrate under the energy distribution curve of Fig. 3. 


From the above relations the dependence of the maximum permissible cur- 
rent for resonance acceleration on the radio frequency power and the Q of the 
cavity may be determined. 


_ 2W, 


| ie 


a ies, (33) 
Vi 

For constant power input, it is of interest to see how the maximum current 
depends on the parameters of the resonator. At constant frequency it may 
be shown that b 


= ’ 

dQ 
where d is the width of the accelerating gap, assumed inversely proportional 
to the equivalent circuit capacity of the resonator. Thus to obtain large 


currents a high Q resonator is desirable. The accelerating gap width, however, 
is primarily determined by other considerations. 


(34) 





oc 
max 


From Equation (33) we find that the maximum orbit current permissible in 
the present accelerator is approximately 0.95 amp., taking W, = 300 kw., 
V = 500 kv., and-G,, as 5 X 1077 mho. 


From an integration under the curve of current against energy the total 
current at all energies in the first orbit is found to be related to J,,,, as 


| = 2.1 Dean 
= 2 amp. 


The effect of the current in the outer orbits is neglected since it is small com- 
pared to the first orbit current. 


The mean current emerging from the cavity in the first orbit is related to 


6 


—) 


. : . 

tmax loons /2 2 
where 6 is the duty cycle = 3.45 K 1074. It will be observed that the duty 
cycle quoted is half the magnetron duty cycle (2 us at 345 p.p.s.).. From oscil- 
lographic observations it was found that current in the orbits starts only when 
the radio-frequency in the cavity has built up to its peak value. The build-up 
time of the cavity is about | uw sec. and thus the current in the orbits exists for 
only 1 wsec. vy is the fraction of the current within the cavity which escapes 
through the accelerating hole. This factor is estimated as 0.8. 


I’ max in the following way 


The loading current within the cavity is in opposite directions during each 
half of the radio-frequency cycle. Since the current leaving the cavity in only 
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one direction is being measured, it is necessary to introduce the factor 2 in 
the denominator. 


Substituting these values we obtain an approximation for the maximum 
permissible electron current 
tmax <= 190 pa. 


Currents greater than 100 wa have been measured emerging from the hole 
in the cavity under conditions similar to those assumed above. Thus it is 
possible to obtain electron currents from field emission alone which are of the 
same order as the maximum permissible currents. This is corroborated by the 
measurement of current obtainable when a tungsten filament was placed at 
one hole in the cavity. In this case when the filament was hot enough to 
produce appreciable emission the outer orbits disappeared, the cavity voltage 
being decreased below the resonance value by the loading of the electron beam. 
Very large currents were obtainable in the first orbit. 

The emission process appears to be a combination of field, secondary and 
photo emission. These currents are largely dependent on the state of the 
surface of the cavity. It has been found that oxidation of the poles of the 
cavity causes a very great increase in emission, whereas almost zero emission 
results from a highly polished and well cleaned surface. 


In an attempt to increase the emission without unduly lowering the Q of the 
cavity, a thin coating of thorium oxide was sprayed onto the cavities in an 
annular ring around the accelerating holes. Copious emission was obtained 
with a bright blue white fluorescence of the thorium oxide. The current was 
unstable and in all cases exceeded the maximum permissible value for re- 
sonance. 

8. The Accelerating Cavity 


The choice of a suitable cavity resonator was influenced by the following 
considerations: 


1. A suitable shape that would have a short accelerating gap, and which 
would not obstruct the small radius orbits. 

2. The cavity must have as high a shunt resistance as possible. 

3. To reduce the problem of sparking at high powers, all surfaces inside the 
cavity must be well rounded. 


An ellipsoid hyperboloid cavity shape was chosen since it fulfilled the above 
conditions, and also because the properties of such a cavity have been calcu- 
lated (1). The final shape of some cavities was modified considerably from the 
original ellipsoid hyperboloid for mechanical reasons. Fig. 5 is a cross-section 
view of one type of cavity; the block is made from solid copper; power is fed 
‘into the cavity through the circular hole (B). The poles of the cavity consist 
of two conical copper spinnings (C) which are clamped into position by the 
plates (D). A knife-edged contact ensures a good electrical connection between 
the cones and the body of the cavity. 
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A cavity turned from a solid block of copper in two halves and the surfaces 
of the central split lapped to an optical finish gave a high Q when the two 
halves were tightly clamped together. Early cavities were machined from brass 
and the inner surfaces silver plated; these cavities did not give as high a Q as 


the copper cavities. 
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Fic. 5. Clamped-cone type cavity. 


A third type of cavity constructed from two copper spinnings brazed together 
around the equatorial plane was also used. These cavities have given the 
highest currents and have all shown measured Q’s of about 9000. 


Since it is extremely difficult to measure the shunt resistance of cavities with 
large orifices, it was decided to measure Q only, and obtain an estimate of the 
shunt resistance from the proportionality between shunt resistance and Q, for 
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Fic. 6. Variation of Q and resonant wave length with the diameter of the orifice in the cavity. 
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cavities of the same shape. A practical value of the shunt resistance was 
assumed to be the same fraction of the theoretical shunt resistance that the 
measured Q was of the theoretical Q. 


A special technique was developed for the measurement of these high values 
of Q (4), and consisted of a circuit for comparing the frequency response curve 
of a tuned radio-frequency circuit of known and variable Q with the frequency 
response curve of the cavity. When the two response curves were equalized, 
the unknown Q was given in terms of the known Q, and the ratio of the two 
frequencies. These response curves were displayed on a cathode-ray tube 
and the instrument was direct-reading in Q. 


Measurements were made of the effect on Q of the size of the holes in the 
cavity through which the orbits pass (see Fig. 6). From this curve it can be 
seen that holes smaller than 0.350 in. diameter have very little effect on Q or 
resonant frequency. A hole diameter of 3/8 in. was used in the final cavities. 


Maximum Q’s of 10,000 were obtained, yielding an estimated shunt resis- 
tance of approximately 2 megohms. Measured Q’s were about 80% of the 
theoretical Q. 

9. The Radio-frequency System 
A diagram of the radio-frequency system is shown in Fig. 7. The tunable 


magnetron (RK-5586) feeds into a rectangular wave guide containing a phase 
shifter for varying the length of line between the magnetron and the cavity: 
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Fic. 7. Wave guide system. 
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this phase shifter consists of a tapered quartz plate which can be moved across 
the guide by means of a vacuum-tight bellows assembly. A side arm on the 
broad side of the guide introduces a variable load which is effectively in series 
with the cavity. Coupling to the cavity is through a circular hole; the diam- 
eter of this hole may be varied by the insertion of metal rings to adjust the im- 
pedance of the cavity presented to the wave guide. 


The equivalent circuit of the radio-frequency system is shown in Fig. 8. 
The effect of a series resistance in stabilizing the operation of a magnetron 
operating into a resonant load has been adequately treated by many authors 
(6), and will not be repeated here. Suffice it to say that a resistive load in 
series with the cavity stabilizes the frequency of operation of the magnetron. 
In our arrangement 50% of the available power is normally fed into the cavity, 
and 50% into the sand load on the side arm of the wave guide Tee. 


LINE 
STABILIZING 
R S RESISTANCE 


MAGNE TRON 


CAVITY 


Fic. 8. Equivalent circuit of radio-frequency system. 


Two properly placed tuning screws in the side arm allow independent ad- 
justment of the resistive and inductive components of the impedance of the 
side arm. ‘The resistive screw allows adjustment of the fraction of the avail- 
able power that is fed to the cavity. The inductive screw is used to cancel 
any small reactive components in the resistive load. A plot on the Smith chart 
of series impedance versus setting of the two controls is prepared from bench 
measurements and used as a guide in adjusting the complete system. 


The wave guide system is evacuated to reduce the problem of sparking, and 
all the controls are brought out through vacuum-tight bellows so that all pos- 
sible circuit adjustments can be made under vacuum. Operation at pressures 
below 3 X 1075 mm. of mercury has been found to be quite satisfactory. 


10. Results and Observation 


The orbits which appear when the cavity is energized and the magnetic field 
adjusted to the resonant value are observed in a movable fluorescent target 
(see Fig. 1). This target is also fitted with a Faraday collector behind a 4 mm. 
wide slit to allow measurement of the electron current in the various orbits. 
Fig. 9 shows the current to the collector plotted against the distance of the 
collector slit from the center of the cavity. The first peak on the right-hand 
side of the diagram is the primary spectrum of energies of the electrons after 
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one acceleration, i.e., the spectrum shown in more detail in Fig. 3. The further 
peaks are the orbits formed after successive accelerations. It will be observed 
that there is little loss of current once the electrons are trapped into the second 
orbit. The variation of the magnetic field in the median plane is shown in 
Fig. 10 and the position of the orbits and cavity are indicated. The accelerating 
gap of the cavity was placed at a point where the magnetic field was 2% lower 
than its value at the center. The eighth orbit then appears at a point where 
the magnetic field is 10% lower than the central value and the rapid decrease 
of the field at this point causes the broadening of the orbit. No shimming of 
the magnet was found necessary to produce adequate focusing of the orbits. 
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Fic. 9. Spatial distribution of current in the orbits. 


The results shown in Fig. 9 are obtained when the magnetic field is adjusted 
to give maximum current in the outer orbits and under the conditions shown, 
i.e., magnetic field = 1000 oersteds the currents are about 1/2 wa. However, 
resonance acceleration with smaller current outputs is observed to occur over 
a range of magnetic field settings, and normally eight orbits are obtained with 
the magnetic field set anywhere between 900 and i200 oersteds. At 1200 
oersteds the energy of the electrons in the eighth orbit is 5.0 Mev. but the 
current is only about 1% of the value at 1000 oersteds. 


Of the electrons leaving the gap after the first acceleration only a small 
fraction have the correct energy and phase to return to the accelerating gap 
in the correct phase for further resonance acceleration. The total electron 
current leaving the accelerating gap is about 100 wa. when the magnetic field 
is set to the resonance value. Of this current, about 1 ya. is trapped into the 
second orbit and subjected to further acceleration. Thus, of the electrons 
originally accelerated only a small fraction (about 1%) are subject to resonance 
acceleration. Since the fraction of the first orbit current which is further 
accelerated is so small, most of the electronic loading of the cavity is caused by 
electrons which are not subject to resonance acceleration. It must be remem- 
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bered that the currents referred to’ above are mean currents and that the 
electronic loading is a function of peak currents which are of the order of 
fractions of an ampere. 
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Fic. 10. Variation of magnetic field within the vacuum chamber. 


Some measurements have been made of the X-radiation produced by the 
4 Mev. and 4.6 Mev. electrons in the seventh orbit and eighth orbit respec- 
tively. For example, when a thick gold target is set to intercept the seventh 
orbit the X-ray dosage in the forward direction measured with a Victoreen 
r-meter 15 cm. from the target is about 20 r. per min. when the orbit current 
is 1/2 wa. 
11. Discussion of Results 


The construction of this accelerator has proved the validity of Veksler’s 
original suggestion, and allows us to consider the advantages and disadvantages 
of this type of electron accelerator. The advantages of this type of accelerator 
might be classified as: 

(1) Extreme simplicity of the accelerator and its associated radio-frequency 

system. 

(2) Production of reasonably large orbit currents without the use of any 
electron injection system. 

(3) Radiation losses from the accelerated electrons at very high energies are 
negligible compared to the energy gained per revolution. Thus the high 
energy limit should far exceed the limiting energy for the betatron. 

(4) Extraction of the beam is a simple matter since the orbits are widely 
spaced from one another. 
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The attendant disadvantages are also enumerated: 


(1) A solid magnet is required, so that the cost of the machine will vary 


ree 


approximately as the cube of the final energy. 

At high accelerating energies, within the limits of present generators, 
and higher magnetic fields the radius of the first orbit will become so 
small that the cavity will obstruct the orbit. Experiments with cavities 
suitably slotted to allow the passage of the orbit have been successful 
and it appears that this problem may be overcome with this type of 
cavity. 
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